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Introduction 


This  progress  report  gives  some  preliminary  observations  of 
the  exploding  drop  phenomenon  which  sometimes  occurs  during' 
burning  of  emulsified  fuels*  It  contains  what  may  be  considered 
a  first  attempt  at  explaining  the  so-called  ^micro- explosion* 
effect  which  has  been  observed  in  the  burning  of  such  fuels. 

The  purpose  of  this  work  was  twofoldi  1)  to  offer  an  explanation 
for  the  droplet  explosion  effect  and  thereby  predict  the  temperature 
at  which  an  emulsified  fuel  drop  would  explode*  and  2)  to  report 
the  results  of  an  experimental  attempt  to  measure  this  maximum 
or  explosion  temperature. 

The  model  for  the  explosion  phenomenon  is  based  on  the  assumption 
that  boiling  at  the  internal  phase-fuel  interface  in  the  emulsion 
is  responsible  for  droplet  break-up.  An  expression  for  the  rate 
at  which  bubbles  form  homogeneously  at  the  interface  is  presented 
and  used  to  predict  the  corresponding  nucleation  temperature » 

A 

This  expression  is  a  very  simplified  form  of  a  description  of 
nucleation  in  systems  of  more  than  one  active  species*  The  assumption 
is  basically  that  the  flow  of  nuclei  past  the  critical  size  is 
independent  of  the  molecular  species  causing  the  transition. 

This  assumption  is  probably  valid  only  for  nuclei  which  are  already 
so  near  to  the  critical  size  that  they  have  substantially  reached 
their  equilibrium  composition.  A  more  accurate  treatment  would  have 
to  take  into  account  the  shape  of  the  free  energy  surface  to  predict 
the  rate  of  flow  of  nuclei  through  the  corresponding  saddle  point. 

The  results  for  these  two  cases  should  not  be  expected  to  differ  too 
greatly  because  the  energy  of  forming  the  critical  nucleus  is  the 


same  for  both  cases. 

The  experimental  method  consisted  of  injecting  a  drop  of  the  test 
emulsion  in  the  bottom  of  a  column  filled  with  a  heavier  immiscible 
liquid  under  a  temperature  gradient.  The  rising  droplet  is  progressive 
heated  until  a  temperature  is  reached  at  which  point  the  droplet 
suddenly  explodes.  This  configuration  models  combustion  of  a 
single  fuel  drop  in  that  the  droplet  surface  temperature  from 
the  point  of  ignition  to  steady  state  burning  increases.  The 
steady  state  droplet  burning  temperature  is  limited  by  the  boiling 
point  of  the  fuel.  However,  in  the  experiment  described  here  the 
emulsion  drop  can  be  heated  to  a  temperature  as  high  as  90#  of  the 
critical  temperature  of  the  fuel.  The  advantages  of  this  are  obvious 
when  one  considers  that  for  some  fuel  emulsions  the  predicted 
nucleation  temperature  is  higher  than  the  boiling  point  of  the  fuel. 

Pour  water-pure  fuel  emulsions  were  tested*  1)  water-decane j  2) 
water-dodecane s  3)  water- te tradecane »  and  4)  water-hexadecane . 

The  predicted  and  measured  nucleation  temperatures  are  in  remark¬ 
ably  good  agreement  (that  is,  with  the  highest  measured  temperatures). 
The  results  indicate  that  for  free  droplet  burning  at  1  atmosphere 
pressure  only  the  water-hexadecane  emulsion  should  burn  ex¬ 
plosively. 

It  must  be  finally  pointed  out  that  it  is  implicitely  assumed 
in  this  work  that  an  emulsion  can  indeed  exist  in  the  superheated 
state  as  a  fuel  drop  containing  relatively  immobile  water  micro¬ 
droplets,  and  that  these  water  drops  are  in  fact  immiscible  in  the 
fuel  even  at  high  temperatures.  This  assumption  has  not  yet  been 
rigorously  tested,  but  it  is  so  basic  to  any  work  concerning  the 


burning  of  emulsified  fuels  (by  definition)  that  it  is  for  the 
present  accepted. 

The  period  of  this  work  was  approximately  from  September  1975 
to  July  1976  and  formed  part  of  a  program  generally  concerned 

1 

with  the  burning  characteristics  of  emulsified  fuels  in  the 
Guggenheim  Laboratories  at  Princeton  University. 
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1.  Nucleation  Rate  and  Equilibrium  Number  of  Nuclei 

The  spontaneous  formation  of  tiny  vapor  nuclei  in  bulk 

liquid  is  a  random  event.  Density  fluctuations  in  the  liquid 

* 

create  cavities  or  voids  which  act  as  sites  for  the  vaporiza¬ 
tion  of  single  molecules.  These  molecules  randomly/  collide  with 
or  escape  from  a  vapor  cavity  and  increase  or  decrease  its 
size  by  a  proportional  amount.  Growth  occurs  by  single  vapor 
molecules  being  absorbed  on  the  surface  of  the  bubble  one  at  a 
time.  The  probability  of  simultaneous  collisions  of  two  or  more 
molecules  is  very  low  and  not  considered. 

Growth  continues  in  this  way  until  a  certain  size,  known  as 
the  critical  size,  is  reached.  Further  increasing  the  number 
of  molecules  in  the  bubble  beyond  this  size  leads  to  a  rapid 
growth  to  macroscopic  size,  often  observed  as  explosive  vapori¬ 
zation.  The  critical  size  represents  the  dividing  point  between 
a  stable  liquid  and  a  stable  vapor  phase.  Right  at  the  critical 
size  the  bubble  is  in  a  metastable  state. 

The  kinetics  of  bubble  growth  up  to  the  critical  size  is 
assumed  to  occur  at  a  slow  enough  rate  so  that  the  principle  of 
detailed  balancing  may  be  employed.  Beyond  the  critical  size 
growth  is  rapid  and  explosive  boiling  occurs.  Growth  beyond 
the  critical  size  cannot  therefore  be  described  by  this  model. 
But  this  is  perfectly  acceptable  because  all  bubbles  which  reach 
the  critical  size  are  considered  to  grow  spontaneously  to  macro¬ 
scopic  size.  2  e  problem  therefore  ends  at  the  point  where  the 


theory  does  not  apply. 


It  is  sufficient  to  have  just  one  nucleus  in  a  given  volume 
reach  the  critical  size  for  explosive  boiling  to  occur.  When 
this  happens  in  an  emulsified  fuel,  the  primary  droplet  may  be 
completely  shattered  by  the  rapid  growth  of  the  nucleus  to  macro- 

i 

sco pic  size.  In  the  general  case  the  number  of  nuclei  which 

grow  to  critical  size  is  usually  termed  the  nucleation  rate.  J, 

in  units  of  nuclei  formed/cm^-sec.  The  nucleation  rate  depends 

very  strongly  on  temperature,  as  expected,  and  is  negligible  for 

all  but  a  very  narrow  temperature  range.  This  range  determines 

the  temperature  of  explosive  boiling. 

These  ideas  have  been  utilized  extensively,  in  the .literature 

to  determine  the  nucleation  rate  for  boiling  in  the  bulk  of  a 
liquid  -  homogeneous  nucleation.  Nucleation  at  a  liquid-liquid 

interface  requires  only  a  physically  based  extension  of  current 

theories. 

Consider  fig.  1.  A  vapor  bubble  has  formed  at  the  interface 

between  liquids  a  and  b.  Unlike  homogeneous  nucleation,  as  for 

example  in  the  bulk  of  liquid  a,  both  a  and  b  molecules  have 

access  to  the  nucleus.  The  nucleus  can  grow  or  decay  by  molecular 

interchange  with  either  a  or  b  molecules.  It  .is  assumed  that 

molecules  of  a  cannot  penetrate  the  interface  and  enter  liquid  b 

(vice  versa  for  b  molecules).  This  is  not  strictly  correct 

because  the  interface  is  really  a  region  in  which  the  liquids  are 

mutually  saturated.  Nevertheless  this  so-called  "rigid  plane" 
model  of  the  interface  between  two  liquids  leads  to  results  which 

are  consistent  with  experimental  measurements  and  so  will  be 

adopted  here. 


With  this  understanding  of  the  interface  it  is  seen  that 
b  molecules  can  only  enter  the  nucleus  by  colliding  with  that 
part  of  the  nucleus  surface  lying  on  the  b  side  (the  heavy  dark 
surface  shown  in  fig.  1),  and  a  molecules  can  only  be  absorbed 
into  the  nucleus  by  collisions  with  the  surface  of  the  nucleus 
lying  on  the  a  side.  These  surfaces  are  called  the  transfer 
area  for  b  molecules,  S^,  and  the  transfer  area  for  a  molecules. 

S  .  respectively. 

Performing  an  analysis  similar  to  that  presented  in  refer¬ 
ence  3  but  now  accounting  for  bubble  growth  by  molecular  ab¬ 
sorption  on  more  than  one  transfer  area  (there  is  only  one  transfer 
area  in  homogeneous  nucleation) ,  the  rate  of  formation  of  nuclei 
per  unit  volume  and  time  which  grow  to  the  critical  size  is 

j  =  - —  1 

Jdx 

.~n&7 
»  - 

where  the  summation  extends  over  molecules  of  kind  i  which  have 
access  to  the  nucleus  through  transfer  area  S^. 

0^  is  the  rate  at  which  molecules  of  kind  i  strike  and  are 
absorbed  on  transfer  area  S^(x).  For  spherical  and  perfectly 
elastic  ideal  gas  molecules  the  number  of  i  molecule;  collisions 
per  unit  of  transfer  area  and  time  is 


^TTm^kT* 

m^  is  the  molecular  mass,  P£  is  the  partial  pressure  of  species  i 


# 


at  the  external  liquid  pressure,  k  is  the  Boltzmann  constant, 
and  T  is  temperature.  Equation  2  accounts  for  the  basic  growth 
mechanism  in  nucleation  theory  -  random  molecular  collisions  on 
the  surface  of  the  nucleus. 

n(x)  is  the  number  of  nuclei  containing  x  molecules  in  a 
unit  volume  when  J  =  0  where 

x  ®  E  x^  3» 

It  may  be  assumed  outright  as  a  Boltzmann  distribution  (ref.  4). 
Alternatively  the  system  of  nuclei  and  liquid  molecules  can  be 
modelled  as  an  ideal  dilute  -solution  and  the  free  energy  of  mixing 
minimized.  In  either  case  the  result  is  the  same* 

n(x)  =  N  exp(-W(x)AT)  4 

where  W(x)  is  the  thermodynamical  work  of  forming  a  vapor  nucleus 
containing  x  molecules,  and  N  is  the  number  of  liquid  molecules 
in  a  unit  volume  with  which  the  n(x)  nuclei  are  "mixed"  to  form 
this  hypothetical  solution. 

The  limits  of  integration  in  equation  1  extend  over  nuclei 
containing  from  1  to  an  infinite  number  of  vapor  molecules  (a 
single  vapor  molecule  is  viewed  as  the  smallest  nucleus).  This 
boundary  condition  arises  from  the  assumption  that  there  are  no 
nuclei  containing  x>>xcri-t  molecules.  It  is  mathematically  con¬ 
venient  to  take  this  number  as  infinity.  This  and  other  approxi¬ 
mations  used  in  homogeneous  nucleation  theory  are  discussed  in 
reference  2.  It  is  assumed  here  that  Uiese  same  approximations 
apply  to  nucleation  at  a  liquid-liquid  interface. 


The  value  of  N  in  equation  4  depends  on  the  location  of 
bubble  formation.  For  nucleation  inside  a  liquid  drop  of  radius 
r  containing  liquid  a»  the  number  of  molecules  in  the  drop  is 

i 

Na  =  <r/ra>3  5a 


where  r  is  the  radius  of  the  assumed  hard  sphere  molecule, 
a 

Eq.  5a  assumes  that  the  total  droplet  volume  is  the  sum  of  e 
hard  sphere  molecular  volumes  of  the  molecules  inside  the  .  uid 
drop.  Similarly  the  number  density  of  a  molecules  on  • 

droplet  surface 

Nas  =  (r/ra>2  5b- 

Equation  5b  assumes  that  the  total  surface  area  of  a  liquid  drop 
of  radius  r  is  tie  sum  of  the  surface  areas  of  the  hard  sphere 
molecules  of  radius  r  ,  Combining  equations  5a  and  5b  *  the 
number  of  a  molecules  on  the  surface  of  the  drop  is 


6. 


The  interface  between  two  liquids  really  consists  of  both 
a  and  b  molecules.  Fig,  2  shows  a  hypothetical  illustration  of 
such  an  interface.  If  the  A-C_region  in  fig.  2  defines  the  inter 
face,  the  number  of  molecules  in  a  unit  volume  at  the  interface 
will  be  the  sum  of  the  number  densities  of  the  a  and  b  molecules 
so  that 


N  *  N 


2/3 

a 


7. 


(The  number  density  of  b  molecules  is  referred  to  a  droplet  of 


-6- 


ft* 


ft* 


ft* 


2/3 

radius  r  containing  only  b  molecules  so  that  ) . 

The  equilibrium  number  density  for  homogeneous  nucleation 

within  liquid  a  can  now  be  written  as  follows  1 
1 

n(x)  =  N  exp(-W(x)/kT)  8. 

SL 


For  nucleation  at  the  interface  between  liquids  a  and  b,  equation 
7  substituted  into  equation  4  gives 

n(x)  =  (N2//^  +  N2//^)  exp(-W(x)/kT)  9. 

The  work  terms  in  equations  8  and  9  are  not  the  same  because 
the  geometry  of  bubble  growth  in  the  bulk  of  a  liquid  is  different 
from  that  at  the  interface.  The  difference  is  primarily  that 
of  a  spherical  nucleus  (homogeneous  nucleation)  and  a  lens 
shaped  bubble  (nucleation  at  the  interface  between  the  two 
liquids) . 

For  nucleation  in  bulk  liquid,  N2^  -*>  N  as  one  must  now 
consider  the  hypothetical  nuclei-liquid  molecule  solution  as 
made  up  of  the  bulk  liquid  molecules  and  vapor  nuclei  with  each 
liquid  molecule  within  the  drop  being  viewed  as  a  potential 
bubble.  This  modification  is  further  discussed  in  reference  5* 

For  nucleation  at  the  interface  between  two  liquids  equation 
1  becomes 


[Sa(x)Sa+Sb(x)Sb] 


1 


00 

dx 

nTxT 


ft* 


10 


This  equation  has  the  form  J  *  Jft  +  where  Ja  is  the  number 

of  nuclei  formed  as  a  result  of  growth  by  absorption  of  only 

a  molecules  on  transfer  area  S  (J^  is  similarly  defined). 

“  a 

Combinirlg  equations  2,4 ,  and  10  gives  the  nucleation  rate 

at  the  interface  between  two  liquids* 

N  xrflfa <f>  +  !i!b 

J  -  TSrST)2  L  m* 

.a  d 


where 

N  =  N2/3  +  N^3 

a  b 

and 

1  ■  *a  + 

2 ,  Specific  Location  of  Bubble  Growth  at  the  Interface 

To  determine  the  nucleation  rate  from  equation  11  it  is 
necessary  to  know  the  work  of  formation  of  a  vapor  nucleus 
containing  x  total  molecules.  This  work  should  obviously  depend 
on  the  geometry  and  particular  location  of  bubble  growth.  For 
a  droplet  of  liquid  a  in  .  .  -o  „  '  liquid  b  a  vapor  nucleus 

can  form  at  any  of  the  following  positions  (see  fig.  3)* 


f; 


dx 


exp(-W(x)/kT) 


3. 


1,  completely  within  liquid  a  (homogeneous  nucleation 
in  liquid  a  -  fig.  3a) I 

2,  at  the  interface  between  liquids  a  and  b  but  entirely 
surrounded  by  liquid  a  except  possibly  for  an 
opening  large  enough  to  permit  b  molecules  to  enter  - 


3* 


between  liquids  a  and  b  and  actually  separating  the 
interface.  The  bubble  can  be  considered  as  two 
spherical  caps  joined  together  -  a  lens  (fig.  3c) » 

4,  at  the  interface  and  entirely  surrounded  by  liquid 
b  except  for  an  opening  just  large  enough  to  permit 
a  molecules  to  enter  the  bubble  (fig.  3d) l  and 

5.  completely  within  liquid  b  (homogeneous 

nucleation  in  liquid  b  -  fig.  3e). 

When  liquid  a  is  more  volatile  than  liquid  b  it  is  necessary 
to  introduce  the  concept  of  a  bubble  completely  surrounded  by 
liquid  b  except  for  an  opening  just  large  enough  to  permit 
a  molecules  to  enter.  Without  this  artifice  a  bubble  forming  in 
this  way  will  contain  only  b  molecules  so  that  the  energy  of 
bufcble  formation  will  be  the  same  as  for  homogenous  nucleation  in 
liquid  b.  If  the  temperature  in  liquid  b  is  unable  to  reach  its 
corresponding  superheat  temperature,  explosive  vaporization  will 
not  occur.  This  will  certainly  be  the  case  in  emulsified  fuel 
combustion  because  the  steady-state  burning,  temperature  does  not 
even  reach  the  fuel  (liquid  b  )  boiling  point  let  alone  become 
superheated.  Yet  in  the  experiments  to  be  reported  the  bubble 
growth  geometry  was  predicted  to  be  the  same  as  that  shown  in 
fig*  3d  and  nucleation  was  still  observed  to  occur  even  though 
the  corresponding  temperature  was  far  below  the  limit  of  superheat 
of  the  fuel.  This  apparent  anomaly  can  best  be  explained  by 
assuming  that  both  a  and  b  molecules  are  in  the  bubble  so  that  the 
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bubble  gas  pressure  is  the  sum  of  the  partial  pressures  of 

a  and  bi  P"  =  P"  +  P".  It  will  later  be  shown  that  the  solution 
—  —  a  b 

of  equation  11  can  be  put  in  the  form 

i 

T  ~  <t3/(P"  “  PQ)2 

where  P  is  the  ambient  liquid  pressure  and  o'  is  surface  tension, 
o 

Prom  this  approximate  proportionality  it  is  seen  that  an  increase 
in  the  gas  pressure  lowers  the  expected  temperature  of  nucleation 
at  constant  PQ  and  c r.  In  the  experiments  which  formed  a  part  of 
this  program,  the  a  and  b  liquids  were  such  that  P£  »  P£  so  that 
the  increase  in  P"  on  including  molecules  of  a  in  the  nucleus 
was  substantial.  This  increase  in  total  gas  pressure  is  apparently 
just  enough  to  bring  theory  and  experiment  into  good  agreement. 

Some  means  must  be  provided  to  enable  molecules  of  a  to  enter 
a  bubble  formed  "entirely"  in  liquid  b,  and  the  concept  of  a 
transfer  area  of  molecular  dimensions  is  introduced.  No  modification 
of  the  theory  leading  to  eq.  11  is  necessary  to  account  for  a 
molecules  in  the  nucleus.  SQ(x)  just  becomes  proportional  to 
the  molecular  surface  area  of  an  a  molecule  , 

When  a  nucleus  is  predicted  to  form  in  liquid  b  away  from 
the  interface  (fig.  3e)  and  P"  »  P£  »  molecular  diffusion  of 
a  molecules  through  liquid  b  must  be  included  in  the  theory 
which  is  based  on  bubble  formation  due  to  random  molecular  motion. 

The  only  thing  which  would  be  effected  if  this  were  done  would  be 
the  pre-exponential  term  in  the  solution  of  eq.  11  , 

These  arguemnts  lead  to  the  conclusion  that  when  one  liquid  is 


v# 


much  more  vo  atile  than  the  other  and  homogeneous  nucleation 
in  the  less  volatile  liquid  is  predicted  to  occur*  the  correspon¬ 
ding  temperature  is  closer  to  what  would  he  expected  for  nucleus 
formation  at  the  interface  -  fig.  3d.  Exactly  how  much  more 
volatile  one  liquid  must  be  than  the  other  for  this  to  be  true 
is  4..  v  unknown  at  this  time.  The  experimental  results 
reported  here  suggest  that  a  one  order  of  magnitude  difference 
in  vapor  pressure  satisfies  this  condition. 

When  a  bubble  forms  entirely  in  liquid  a  as  shown  in  fig. 

3b*  the  above  arguments  apply  if  P£  »  P".  Otherwise  if  P"  »  P 
they  are  superfluous.  Even  if  molecules  of  b  are  in  the  bubble* 
the  much  greater  volatility  of  liquid  a  implies  that  the  gas 
pressure  is  almost  entirely  due  to  the  partial  pressure  of  a. 

The  only  difference  between  bubble  formation  as  shown  in  figs, 

3a  and  3b  is  that  one  must  now  consider  liquid  molecules  at  the 
interface  and  not  liquid  molecules  in  the  bulk  of  liquid  a 
(fig.  3a)  as  potential  nuclei.  The  pre -exponential  factor  in 
the  expression  for  the  nucleation  rate  is  thus  reduced  by  about 
106  (N-*N2//3). 

It  is  now  apparent  that  one  needs  to  know  both  the  location 
and  work  of  formation  of  a^nucleus  to  be  able  to  solve  eq.  11 
because  W(x)  will.be  slightly  different  for  the. various  positions 
of  bubble  growth  shown  in  fig.  3.  For  this  purpose  it  would  be 
useful  to  prescribe  conditions  which  could  be  used  to  estimate 
the  position  of  bubble  formation.  Such  conditions,  if  based  on 
the  physical  properties  of  the  two  liquids,  would  provide  a 


O'  a 


powerful  prediction  tool  for  nucleation  at  a  liquid-liquid 
interface.  These  criteria  can  be  determined  by  first  considering 
the  energy  to  form  a  bubble  separating  the  interface  between 
liquids  a  ancl  b  (fig.  3c).  From  this  geometry  the  limiting 
case  of  a  completely  spherical  nucleus  and  corresponding  conditions 
for  its  attainment  can  be  determined. 

II.  Thermodynamics  of  Bubble  Growth  at  a  Liquid-Liquid  Interface 

1.  Basic  Assumptions 

Bubble  growth  at  a  liquid-liquid  interface  is  assumed  to 
occur  under  the  following  conditions; 

1.  Bubble  growth  up  to  the  critical  size  occurs 
reversibly  so  that  it  is  at  all  times  in  mechanical 
equilibrium^ 

2.  The  temperature  is  constant; 

3.  The  external  liquid  pressure  is  constant; 

4.  The  total  volume  of  nucleus  and  surrounding  liquid 
is  constant' (closed  system); 

5.  The  characteristic  nucleus  dimension  (e.g. ,  radius) 
is  much  smaller  than  the  radius  of  curvature  of  the 
interface  at  which  the  bubble  grows; 

6.  Bulk  liquid  and  vapor  properties  can  be  used  to 
describe  the  thermodynamic  state  of  the  bubble; 

7.  The  vapor  in  the  bubble  is  an  ideal  solution  of 
ideal  gases; 

8.  The  two  liquids  on  either  side  of  the  interface  are 
completely  immiscible;  and 

9.  One  liquid  is  much  more  volatile  than  the  other  so 
lhat  P"  »  Pi  (a  denotes  the  volatile  liquid). 

The  first  assumption  implies  that  Laplace's  equation  can 
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L. 


be  used  to  relate  the  gas  pressure  to  the  bubble  radius.  This 
has  been  critized  in  reference  6  because  no  account  is  taken  of 
viscosity  and  inertia  forces  during  bubble  growth  to  the  critical 

I 

size.  The  equation  of  motion  of  bubble  growth  should  really 
be  solved  to  determine  the  relation  between  the  bubble  gas  pres¬ 
sure  and  bubble  radius.  This,  however,  greatly  complicates 
the  solution  and  would  in  any  case  have  a  small  effect  on  the 
nucleation  temperature  because  such  non-equilibrium  effects  would 
appear  only  in  the  pre-exponential  term  in  the  solution  of  eq.  11. 

The  second  assumption  arises  from  the  fact  that  the  nucleation 
rate  is  negligible  until  a  temperature  range  is  reached  at  which 
the  rate  becomes  very  large  -  the  temperature  range  of  expected 
nucleation.  This  temperature  range  is  typically  within  1C  to  2C 
of  the  limit  of  superheat  for  homogeneous  nucleation  and  should 
also  be  representative  of  nucleation  at  a  liquid-liquid  interface. 
The  entire  event  of  bubble  growth  to  macroscopic  size  can  therefore 

be  considered  to  begin  and  end  (growth  to  macroscopic  size)  at 
a  temperature  which  is  within  about  2C  of  the  actual  superheat 

temperature . 

A  hypothetical  system  is  constructed  consisting  of  a  single 
bubble  and  surrounding  liquid  bounded  by  a  control  surface 
across  which  no  mass  flows  (fig.  4b),  Bubble  growth  proceeds  by 
vaporization  of  the  liquid  molecules  within  this  closed  system. 

The  total  number  of  molecules  is  therefore  the  same  before  and 
after  the  appearance  of  the  nucleus  (figs.  4a  and  4b). 

In  emulsified  fuels  the  internal  phase  diameter  can  be  as 


small  as  10"^  cm.  However,  because  a  nucleus  of  critical  size 
is  typically  10“^  cm  in  diameter  droplet  curvature  may  therefore 

•mil  .7 

be  neglected  (10~  »  10  '). 

The  use  of  bulk  liquid  properties  to  describe  the  state  of 
a  microscopic  bubble  was  studied  in  reference  7.  It  was  concluded 
that  because  the  effect  of  bubble  curvature  on  surface  tension 
cannot  be  accurately  estimated,  it  is  a  good  first  approximation 
to  use  the  bulk  liquid  surface  tension  in  describing  the  state 
of  the  critical  nucleus.  The  experiments  reported  here  and  in 
the  literature  indicate  that  this  is  an  excellent  assumption. 

The  effect  of  vapor  non-ideality  in  homogeneous  nucleation 
was  examined  in  reference  8.  It  was  concluded  that  for  boiling 
in  pure  liquids  this  effect  is  very  small. 

In  addition  to  being  satisfied  for  the  particular  a  and 
b  liquids  used  here,  the  last  assumption  provides  a  mathematical 
convenience  which  enables  eq.  11  to  be  integrated.  Without 
this  assumption  an  analytical  solution  to  eq.  11  cannot  be  ob¬ 
tained  . 

The  reversible  work  of  forming  a  nucleus  at  a  liquid-liquid 
interface  can  now  be  determined  using  the  above  assumptions. 

For  this  purpose  the  following  section  reviews  those  aspects 
of  thermodynamics  which  are  important  to  this  problem. 

2,  Thermodynamics  of  Nucleation 

The  first  law  of  thermodynamics  for  a  reversible  process  is 


dU  =  dQ, 


12 


where  £  dW  represents  all  forms  of  reversible  work  which 
rev 

the  system  can  exchange  with  its  surroundings ,  dU  is  Ihe  change 

of  accumulated  energy  of  the  system,  and  dQ  represents  the 

exchange  of  heat  btft«i<«nthe  system  and  its  surroundings.  For 

bubble  growth  at  a  liquid-liquid  interface  £  <*W  includes 

rev 

contributions  from  expansion,  surface,  and  chemical  work  so  that 
the  total  work  is 

S  dWrev  "  2  -  2  Vxi  ’  •=  0 

u^  and  x^  are  the  chemical  potential  and  number  of  molecules  of 
species  i  respectively. 

The  first  term  in  eq.  13  represents  the  expansion  work  con¬ 
tribution  to  the  total  work  of  a  system  which  is  divided  into 
k  subsystems  each  of  volume  and  pressure  P^.,  If  for  example 
one  of  these  subsystems  consists  of  a  vapor  volume  containing 
two  species  then  Pfc  =  £  P^  =  P"  +  P£.  Writing  the  expansion  work 
term  in  the  form  shown  in  eq,  13  is  convenient,  though  somewhat 
unusual,  for  studying  the  thermodynamics  of  the  nucleation 
process. 

The  summation  in  the  second  term  in  eq.,  13  extends  over  the 

number  of  different  species  in  each  phase  present  in  the  total 

system.  The  energy  required  to  create  or  destroy  surfaces 

within  the  total  system  is  represented  by  the  last  term  in  eq.  13 

where  tf.  is  the  surface  tension  of  the  A.  surface.  The  summation 
3  J 

extends  over  the  number/. surfaces  in  the  system.  These  surfaces 
include  bubble  surfaces  and  interfacial  areas  between  liquids  (in 


which  case  d.  is  the  interfacial  tension  between  the  two  liquids) 
J 

The  second  law  of  thermodynamics  for  a  reversible  process 


'  dQrev  =  Tds  14  • 

The  combination  of  equations  12,13,  and  14  gives  the  total 
internal  energy  in  each  differential  piece  of  matter i 


dU  =  TdS  -  2  P.  dV.  +  2  u.  dx.  +  2  o.dA.  15. 

K  K  1  1  J  J 

When  the  system  for  which  eq.  15  is  valid  undergoes  a  change 
in  size  only,  the  extensive  properties  S,  V,  x,  and  A  will  change 
by  a  proportional  amount.  The  amount  of  these  changes  is  such 
that  the  intensive  properties  T,  P^,  u^,  and  do  not  chnge. 
Under  Viese  conditions  integration  of  eq.  15  between  the  initial 
and  final  states  of  this  hypothetical  size  change  gives 

/W-TAS-SF^  +  tu^Xj+SOjAAj. 

Only  the  size  of  the  system  has  been  increased  so  that  the 
final  state  must  be  some  multiple  of  the  initial  state.  Since 


,  '  extensive  properties  depend  only  on 


size, they  all  change 


by  the  same  amount.  Therefore  from  the  above  equation, 


U  =  TS  -  2  P.  +  2  u-x.  +2  O’. A.  16. 

<*£  11  J  J 

Equation  16  gives  the  internal  energy  of  a  system  of  finite  size, 
The  differential  of  equation  16  is 

dU  =  TdS  +  SdT  -  2  PkdVk  -  2  V^dP^  +  2  ui<ixi  +  2  xidui 

+  2  bjdAi  +  2  A. do.. 


-J.O- 


Ihis  equation  must  still  satisfy  equation  15.  This  implies  that 
eq.  17  below  must  be  valid  1 

SjdT  -  £  VkdPk  +  £  xj_dui  +  2  A^do^.  =0  17 . 

Eq.  17  is  a  form  of  the  Gibbs-Duhem  equation.  For  bubble  growth 

at  constant  temperature,  dT  =  0  and  do-  »  0  since  cr*  is  a  strong 

J  J 

function  of  temperature  and  may  only  very  weakly  depend  on 
pressure.  Eq.  17  then  becomes 

£  VkdPk  =  £  xidui  -18. 

Assuming  each  liquid  and  vapor  species  is  confined  to  only 
one  sub-system  volume,  the  subscript  "k"  in  eq;  18  .can  be  re¬ 
placed  by  the  subscript  "i"  so  Hi  at 

£  V^dP^  =  £  x^du^. 

If,  for  example,  the  ith  species  is  a  vapor  contained  in  a  volume 
,  then  P^  would  be  its  corresponding  partial  pressure.  Under 
the  condition  stated  above,  this  particular  species  in  the  vapor  phase 
will  not  be  present  in  any  other  sub-system  volume. 

It  is  now  assumed  that  for  the  above  equation  to  be  satisfied 
each  term  on  the  left  hand  side  must  be  equal  to  the  term  cor¬ 
responding  to  the  same  species  and  phase  on  the  right  hand  side. 

If  this  is  a  valid  assumption,  one  can  write 

VidPi  =  xidui  19. 


Equation  19  relates  the  Chemical  potential  of  the  ith  species 
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in  the  liquid  (')  or  vapor  (")  phase  to  the  volume  in  which  it 
is  contained  at  corresponding  pressure  P^.  For  an  ideal  gas, 

(  P£V£  =  x£kT 

so  that  eq.  19  becomes 

kT  f£i  =  du? 

Pi 

Integrating  the  above  equation  between  some  arbitrary  aib-critical 

vapor  pressure,  Pv,  and  the  vapor  pressure  in  the  critical 

size  bubble,  PP  ,  gives 
ic 

u?(P£,T)  -  uPc(P?c,T)  =  kT  ln(PP/PVc)  20a. 

A  nucleus  at  the  critical  size  is  in  chemical  equilibrium 
so  that  the  vapor  chemical  potential  of  species  i  is  equal  to 
its  liquid  chemical  potential i 

u£c<P£c»T>  =  u{(P»,T). 

Combining  the  above  condition  with  eq,  20a  relates  the  chemical 
potential  of  liquid  and  vapor  in  a  nucleus  of  sub-critical 
sizes 

u£(P£,T)  -  u?(P£,T)  =  kT  ln(P£/P£c)  20b. 

Equations  1 6  and  20b  will  be  particularly  useful  in  deter¬ 
mining  the  energy  of  formation  of  a  bubble  at  a  liquid-liquid 
interface.  For  this  purpose  it  is  convenient  to  use  the  Helmholtz 
function  because  the  minimum  work  in  a  constant  temperature  pro- 


cess  is  equal  to  the  decrease  of  the  Helmholtz  function  for  the 
system. 

The  Helmholtz  free  energy  is  defined  as 

P  8  U  ■  TS  21. 

The  differential  of  equation  21  is  the  following i 

dF  =  dU  -  TdS  -  SdT  22. 

Substituting  equations  12  and  14  into  eq.  22  gives, 

dF  -  -  £  dWrey  23. 

The  Helmholtz  free  energy  is  therefore  minus  the  total  work 
obtainable  from  an  isothermal  process. 

Equations  16  and  21  give  the  Helmholtz  free  energy  of  a  system 
of  finite  sizei 

F  ■  - £  pk\  +  £  “ixi +  £  Vi  2"- 

This  is  the  basic  equation  used  to  determine  the  free  energy 
of  bubble  formation.  To  do  this,  two  systems  of  finite  size 
are  considered,  one  with  and  the  other  without  a  nucleus.  The 
change  in  free  energy  between  the  two  systems  as  determined  .by- 
eq.  24  applied  to  each  gives  the  minimum  work  of 

forming  a  bubble  of  specified  size. 

3.  Work  of  Bubble  Formation  at  a  Liquid-Liquid  Interface 

Consider  the  total  system  shown  in  fig,  4a  bounded  by  a 
control  volume  (dotted  line).  The  total  system  consists  of  two 
pure  liquid  sub-systems  which  are  separated  by  an  interface 


(modelled  as  essentially  a  dividing  line  between  the  two  liquids). 
The  total  system  consists  of  x^  molecules  of  liquid  a  at  a 
liquid  pressure  P^,  volume  V^,  and  liquid  chemical  potential 
(similarly  for  liquid  b ) . 

After  the  appearance  of  a  single  bubble  at  the  interface, 
the  total  system  consists  of  three  sub-systems  as  shown  in  fig.  4b. 
The  total  system  boundary  is  the  same  as  in  fig.  4a  so  that  the 
total  volume  has  not  changed.  The  presence  of  the  bubble  implies# 
however,  that  the  sub-system,  volumes  must  have  changed.  The 
temperature  is  the  same  before  and  after  the  appearance  of  the 
bubble. 

The  Helmholtz  free  energy  of  the  system  without  the  bubble 
(fig.  4a)  is  the  following! 

Fa  “  xalual  +  xbl“bl  +  8abSabl  '  pklykl  '  PblVbl  25- 

Similarly  for  the  system  in  fig.  4b  the  free  energy  as  determined 
from  eq.  24  is. 


a2  a2 


+  x*  u’  -P* V*  -  P, 


b2  b2  a2  a2 


b2  b2 


ab  ab2 


+  (*b»b  >  xaua  -  P"Vb  +  ffasa  +  CTbsb> 

The  quantity  in  brackets  in  eq,  26  represents  the  energy  contri¬ 
bution  froro the  bubble  at  the  interface.  The  important  point  to 
note  in  eq.  26  is  that  the  presence  of  both  a  and  b  molecules 
in  the  bubble  has  been  accounted  for  by  including  the  chemical 
energy  ■  of  b  molecules  in  the  bubble,  x£u£.  This  is 

consistent  with  the  assumption  that  bubble  growth  occurs  by  molecular 


absorption  on  the  appropriate  transfer  areas,  S&  or  so  that 
P"  “  +  P£  in  eq.  26. 

Including  terms  such  as  cTaAa  or  <rb  A  b  in  eq.  26, as  for  example 
one  might  be  tempted  to  do  to  determine  the  energy  of  bubble  for¬ 
mation  for  growth  as  shown  in  figs.  3b  and  3d,  will  always  be 
negligible  compared  to  either  cr&Sa  or  obSb  except  possibly  for 
conditions  which  would  be  of  no  practical  interest. 

The  following  relations  apply  as  a  consequence  of  the 
previously  mentioned  assumptions i 

*7l 
28, 


xal  =  xa2  +  xa 


*61  =  xb2  +  xb 


al 


+  Vbl  -  +  VKP)  = 


a2 


b2- 


29. 


The  number  of  a  and  b  vapor  molecules  is  assumed  to  be  small 
compared  to  the  number  of  a  and  b  liquid  molecules  respectively. 
With  constant  a  and  b  liquid  pressures  (assumption  3  on  page  11), 
the  liquid  chemical  potentials  in  the  two  systems  shown  in  figs.  4a 
and  4b  will  be  the  same  j 


u*,  »  u* 


al 


a2 


ubl  w  ub2 


30, 

31. 


In  the  specific  application  of  the  results  of  this  inves¬ 
tigation,  liquid  a  is  dispersed  as  small  drops  in  a  larger  mass 
of  liquid  b.  There  will  be  a  finite  curvature  of  the  interface 
between  a  and  b  (see  fig.  5).  Because  the  drop  of  liquid  a  is 
in  mechanical  equilibrium,  the  pressure  across  the  interface  is 


given  by  the  Laplace  equation i 


Pal-  p£l  -  2oab//rab  32 

where  cr&b  and  tr&b  are  liquid  interfacial  tension  and  liquid 
droplet  radius  respectively.  For  conditions  of  present  interest 
2cr&l3/rat  is  much  less  than  either  orlP'b^.  The  pressure  in  the 
drop  and  surrounding  liquid  are  therefore  approximately  the  same 


(assumption  3) * 


pal  - 


Pa2  w  Pb2* 


For  bubble  growth  at  constant  liquid  pressure  it  follows  that 
Pll  -  Pbl  -  Pk2  -  P^2  =  Po‘ 

Subtracting  eq.  26  from  eq.  25  gives  the  total  system  free 
energy  change  on  bubble  formation i 


xalual  +  ^b^bl  +  crabSabl  "  PalVal 

“Xa2Ua2  "  HzHz  +  Pa2Va2  +  Pb2Vb2 
-  x"u»  +  P"V£  -  <?aSa  -  CbSb 


PblVbl 


aabSab2 


Grouping  terms  and  setting  all  liquid  pressures  equal  to  F  ,  eq.  33 
becomes, 

Fa  "  Fb  "  P"VJ  ’  Po<Yil  +  vil  -  ^a2  "  Vfe  )  -  “  °bSb 

"  °ab^Sab2  "  Sabl}  +  xalual  *  xblubl  “  xa2ua2 
‘  HzKf  “  xbub  -  xaua 


34 


Substituting  equations  27  to  29  into  eq.  34  and  grouping 
terms  gives  the  following 1 

pa  “  Pb  Vb^P"  "  “  °asa  *  °bsb  *  °ab<sab2  '  Sabl^ 

*  xj2<ub2  -  »4l>  *  xLzK2  -  “al>  -  xb<ub  -  ubl- 

-  xa(ua  -  “il1 


With  equations  30  and  31 »  eq.  35  becomes 


Pa  -  Pb  "  Vb(P"  -  V  -  °aSa  *  aab(Sab2  *  W 


-  -  “bl5  -  xa(uS  ‘  ual> 


36. 


The  interfacial  area  between  liquids  a  and  b  destroyed  in 
forming  the  vapor  bubble  (fig.  4b)  is  related  to  S&b2  and  S&bl 


simply  as 


si  =  Sabl  *  Sab2-  Eq.  36  can  then  be  re-written  as 


F  -  F  =  V"(p«  -  p  )  -  a  s  -  c^s,  +  0  .  s. 
a  b  bv  o'  a  a  b  b  ab  i 


-  x£(u£  -  u^)  -  x:(u:  -  ul-,) 


37. 


bl 


a'  a 


al; 


From  eq.  23  the  minimum  work  required  to  form  a  vapor  nucleus  at 
the  interface  between  two  liquids  is 


W  =  x£(u£  -  U£x)  +  x£(u£  -  njp  +  V£(p0  -  P«) 


+  daS_  +  CyS.  -  OaKS. 
a  a  b  b  ab  1 


38. 


The  chemical  potential  differences  in  eq.  38  are  related  to  the 
component  partial  pressures  through  eq.  20b  so  that  eq,  38  can  be 


re-written  as  follows i 


W  *  x£kT  ln(P£/PJc)  +  x£kT  ln(P£/P”c)  +  V£(Po  -  P") 


+  c  S  +  d,  S  -  a  .  S. 
a  a  b  b  ab  x 


39. 


The  equation  of  state  of  an  ideal  gas  can  be  used  to  eliminate 
x£  and  x£  in  eq.  39  (i.e.,  x£  =  PJV^/kT  and  x"  =  P£V£/kT)i 


W  =  VJ(F£ 


lnCPj/P^)  +  P”  ln(P"/P"c) 
+  °aSa  +  ®bSb  ‘  ffabSi 


+  P 

o 


where  P"  in  eq.  40  has  been  replaced  by  P£  +  P£. 

Eq.  40  is  the  final  form  of  the  work  to  form  a  vapor  bub¬ 
ble  at  the  interface  between  two  liquids  under  the  previously 
described  assumptions.  The  first  two  terms  in  parenthesis 
represent  chemical  work,  the  third  term  is  the  work  directed 
against  pressure  forces,  the  fourth  and  fifth  terms  represent 
the  work  to  create  the  vapor  surface  areas  S&  and  s^,  and  the 
last  term  is  the  work  to  destroy  the  liquid-liquid  interfacial 
area  . 

When  liquid  a  is  much  more  volatile  than  liquid  b,  P£  »  P£ 
and  eq.  40  reduces  to  the  expression  given  in  reference  4. 

For  simplicity  it  is  assumed  that  the  nucleus  shown  in  fig.  4b 
is  a  lens  comprising  two  spherical  caps  joined  together. 

With  reference  to  fig.  6a,  the  bubble  volume,  V£,  can  be  ex¬ 
pressed  in  terms  of  r  ,  r.  ,  and  the  contact  angles.©  and  4>  . 

From  geometry  the  volumes  of  the  two  spherical  segments  shown  in 
fig.  6a  are, 


4] 


=  l/3tth;(3ra  -  ha)  =  l/3lTr^(3(ha/r &)c  -  (h/V  } 

and 

vj  =  l/3^r3(3(hb/rb)2  -(h,/rb)3)  41D. 

where  1 

V "  =  v "  +  v"  -Ale. 

b  b  a 

Again  from  geometry, 

M  =  cos(t)=  (r  -  h  )/r  =  1  -  h/r  42a 

a  a  a  a  a  a 

and 

=  cos(a)  =  1  -  h^/r^  42b. 

Combining  equations  41  and  42  and  rearranging  terms,  the  total 
bubble  volume  is 

V£  =  1/3  Tf  (z^(2-3Mb+M^)  +  r^(2-3Ma+M-J))  43. 

For  the  nucleus  to  be  in  mechanical  equilibrium  the  forces 
at  the  edge  of  the  lens  must  balance,  A  force  balance  at  the 
bubble  edge  shown  in  fig.  6b  results  in  the  following  expressions i 

o  ,  -  o'  cos (9)  -  o',  cos  (4>)  =  0  44a j 

ab  a  b 

and 

0  sin(9)  -  c,  sin(4)  *  0  44b. 

a  o 

Combining  equations  44  and  and  using  the  identity 
2  2 

cos  (9)  +  sin  (9)  =  1,  one  can  also  write  the  followingi 


cos(S)  »  (®a+0ab-0b)/2<Vat> 

45a  i 

oosO)  =  (<^+®ab-°a)/2a'b0ab 

45b. 

Because  the  line  of  force  of  surface  tension  is  tangent 


and 


to  the  spherical  segment  at  the  apex  of  the  bubble  and  forms 
a  right  angle  with  the  respective  radii  of  curvature  at  that 
point ,  it  follows  that  9  =  ^  and  a  =  (j>  .  Therefore,  and 

can  be  expressed  in  terms  of  the  surface  and  interfacial 
tensions  by  combining  equations  42  and  45* 

“a  =  cos(6)  =  (aa+sab-®b)/2oa<rab  46a' 

and 

«b  =  C0SW  =  <0b+«a1J-«a)/2Vab  46b- 


These  expressions  were  derived  in  references  4  and  9.  They  hold 
theoretically  only  at  the  critical  size,  but  a  negligible  error 
can  be  expected  by  applying  them  to  nuclei  which  are  "near" 
the  critical  size  as  mentioned  previously. 

The  radii. of  curvature  of  the  shperical  caps,  r&  and  r^, 
are  related  to  each  other  as  follows i 

sin(a)  =  sin(<}>)  =  i/r^j 


and 


sin  (S')  =  sin(9)  =  a/r 

CL 

where  a  is  defined  in  fig.  6a  and  may  be  considered  as  the  charac 
teristic  lens  dimension.  Combining  the  above  two  expressions 
provides  a  formula  relating  a  to  r&  and  r^i 

a  =  r  sin(9)  =  r.  sin(<J>)  47. 

a  b 


Eliminating  sin(9)  and  sin(<}>)  between  equations  47  and  44b  gives 
the  relation  between  radii  of  curvature  and  surface  tension i 

ra/oa  =  Vffb  48  • 


The  bubble  volume  can  now  be  expressed  in  terms  of  either  r  or 


x*  by  substituting  equation  48  into  equation  43 1 

D 


V£  =  ((2-3Ms+m|)  +  (2-3Mb+M^)o^/^)irr^/3  49a* 

or  alternatively 

V£  =  ( (2-3M^+M^)  +  (2-3Ma+M^)c^/<y^)lTr^/3  49b. 


By  defining 

F  =  (^(2-3Ma+M^>  +  ff^(2-3Mb+M^))/4 

and 


P  =  °aFa  =  °bPb 
the  bubble  volume  becomes, 


V"  =  r^F  ff  4/3  =  r^F.  4/3 
D  a  a  d  d 


50a 

50b, 

51. 


F  and  F  in  equations  50  and  51  can  now  be  considered  as 

Cl  D 

corrections  applied  to  a  spherical  bubble  of  radius  r  or  r  so 

a  b 

that  its  volume  equals  that  of  a  lens  whose  surfaces  have  radii 
of  curvature  r&  and  r^.  It  should  therefore  be  expected  that 
for  a  spherical  bubble  of  radius  r  ,  F_  =1  (figs.  3a,  3b,  3d,  and 

Cl  3 

3e),  The  conditions  under  which  F„  =  1  and  F,  =  1  will  be 

a  o 

shortly  described. 

The  surface  areas  of  the  spherical  caps  are, 


and 

where 


s*  =  4l,rbAb 

A.  =  <l-Ma)/2 

^  “  d-v/2 


52a, 

52b 

52c 


and 


52d 


The  liquid -liquid  interfacial  area  destroyed  by  growth 

of  the  vapor  lens  of  "radius"  a  is  =  TT  a  ,  Combining  this 

2  2 

with  equations  46a ,  47,  and  the  identity  cos  (©)  +  sin  (©)  =  1» 


the  interfacial  area  becomes 

Si  "  «  ra(l-«4> 


52e. 


The  work  to  form  a  bubble  at  the  interface  between  two 
liquids  can  now  be  expressed  as  a  function  of  r&,  P£,  and  P"^ 
by  substituting  equations  51  and  52  in  eq.  40  and  rearranging 
terms i 

W  =  (p;  ln(pj/p;0)  +  PJ  ln(P£/PJc)  ♦  P„-  PJ  -*{>*^4/3 

+  rfaaF 
a  a  a 


At  the  critical  size,  P"  =  P"  and  P"  =  P"  so  that  terms 

a  ac  b  be 

containing  the  natural  logarithm  in  eq.  53  will  vanish.  If  the 
bubble  is  near  enough  to  the  critical  size  for  these  terms 
to  be  very  small  compared  to  P  in  eq.  53,  it  should  be 

expected  that  a  Taylor's  series  expansion  of  eq.  53  about  cor - 
ditions  at  the  critical  size  will  be  a  good  approximation  to  the 
work  of  bubble  formation.  In  addition  the  dominant  contribution 
to  the  integral  in  eq.  11  comes  from  conditions  at  the  critical 
size  (the  extension  of  the  limits  of  integration  from  -od  to  +cd 
is  made  because  of  this  fact).  The  accuracy  of  this  and  other 
approximations  used  in  nucleation  theory  is  discussed  in  refer¬ 
ences  2  and  10.  It  should  nevertheless  be  pointed  that  an 
analytical  solution  to  eq,  11  can  crilybe  determined  by  using 
the  Taylor's  series  form  of  W  to  perform  the  required  integration. 


W  is  a  function  of  three  variables  -  r„ ,  P" ,  and  P”  -  so 

a  a  b 

that  the  Taylor's  series  appropriate  to  a  function  of  three 
independent  variables  must  be  used  in  the  expansion.  Expanding 
about  conditions  in  the  critical  size  nucleus,  the  three  variable 
Taylor's  series  is 

*  *  W(rac’Pac’Pbc>  +  (r-r^W/Srl  +  (P^c»W/3Pblc 


+  (SS-P^feW/SP-^  +  (>rac)(PJ-P"c)i^P^|t 
+  (P^P“c)(r-rac)92W/^Ppra|c  +  (P^c)  (Pa“Pac^2w^P^P; 

+  ^(r~rac)2')2w/<}ra|c.  +  ®(P"-P^c)2d2W/aP^2jc  54. 

+  i(P^-P|;c)232W/dP£2|t  +  higher  order  terms 

At  the  critical  size  the  bubble  radii  are  related  to  the  f 
pressure  difference  across  the  bubble  surface  as  follows* 


P?  +  P"  -  P_  =  2o  /r  =  2 a,/r.  * 

be  ac  o  a  ac  b  be 

which  is  the  condition  of  mechanical  stability  at  the  critical 
size  (Laplace's  equation).  Using  eqs.  55  and  53,  the  required 
partial  derivatives  in  eq,  54  are  the  following* 


W|c  =  4/3"rac  fa 
J2W/0r^j=  ”8trtfaFa 
0 

)2W/5P£3PJ|  =  0 

)ViPJ2|=  V3irr3cFa/p;c 


uf/fcp;|-  o 

}*W/ir  >PJ|  -  0 


i2W/iP-2|  =  4/3  If  ra0Fa/PaC 

lw/JPSlt“  0 

52w/JraJP;j=  0 


56 


s  the 


W  =  it 


W  must  be  expressed  as  a  function  of  only  r  to  evaluate  the 

cl 

first  term  in  brackets,  and  a  function  of  only  r^  to  evaluate 
the  second  term  in  brackets. 

The  equation  of  state  of  the  gas  in  the  bubble  is 

X"  =  (P"  +  P")V"/kT. 
a  b  b 

Substituting  equations  51  and  55  into  the  above  equation  gives 


the  f olio wing « 


x-  -  V3tTFa(r^Po+2oar^)/kT 


Differentiating  eq.  59  witn  respect  to  r  and  rearranging,  one 

cL 


obtains 


dx"/dr  =  WF  rf(P"+P”)B/kT 

E  E  a  a  D 


where  £  is  defined  as 


B  »  l-(l-P0/(P"+P^))/3  61. 

Expressing  the  equation  of  state  as  a  function  of  r^»  one 
obtains 

x-  «  4/31TFb  (r^^r^/kT  62. 

Differentiating  eq.  62  with  respect  to  rb  gives, 

dx"/drb  »  4TTFbrb(P«+P")B/kT  63. 

To  perform  the  required  integration  in  eq.  58,  P£  and  P£ 
must  each  be  expressed  as  a  function  of  r  or  r^  so  that  W  in 

a  a 

eq.  57  will  be  a  function  of  only  one  variable  (r  or  rb). 

There  are,  however,  four  variables  and  only  two  independent 
relations  -  equations  48  and  55.  Unless  an  additional  equation 
is  available  relating  these  four  variables,  eq.  58  cannot  be 
evaluated  analytically.  It  is  now  apparent  why  in  a  previous 
treatment  of  this  problem  (ref.  4)  the  assumption  of  only  one 
volatile  liquid  was  made,  for  this  provides  the  additional  con¬ 
dition  needed  to  solve  eq.  58. 

When  P"  »  P£,  eq.  55  can  be  written  as 

P-  -  Po  „  P-  +  p£  -  Po  =  20a/ra  =  2<Vrb 

At  atmospheric  pressure  P”  will  usually  be  much  greater  than 

a 

P  .  However,  if  nucleation  occurs  at  high  pressure,  PQ  and  P 
may  both  be  the  same  order  of  magnitude.  For  this  reason  PQ 
will  be  retained  in  eq,  64  with  the  understanding  that  it  has 
a  small  effect  at  P  =latm,. 


m  * 


Near  the  critical  size  where  P”  and  P"c  are  approximately 
the  same  order  of  magnitude  (similarly  for  P£  and  P£c),  "the  con¬ 
dition  P"  »  P"  implies  that 
a  o 


P"2/P" 
a  '  ac 


» 


p-Vp 


be* 


The  second  term  in  brackets  in  eq.  57  can  therefore  be  neglected 
compared  to  the  first.  With  this  simplifying  assumption  eq.  57 
can  be  combined  .with  equations  48,  50,  and  64  to  express  W  in 
terms  of  the  single  variable  r&  or  r^. 

Since , 


2/3tr'Vac(Pa-PSc)2/PM  =  V3ffFaaa(ra-rac)2U-Po/p-c)  65. 


eq.  57  can  be  expressed  as  a  function  of  only  r 

a 

follows  I 

W(r  )  =  4/3 ffr2  o’  F  -  (r  -r  )24tTo  P  B 
a  ac  a  a  a  ac'  a  a  a 

and 

*(*b>  =  4/3ffrbeVb  *  'VV^bVa 

In  equations  66  and  67  B&  «  B  because  P^  »  P£.  The  subscript 
"a**  will  therefore  be  dropped  when  re-writing  equations  66  and 
67  in  the  future. 

After  substituting  equations  52a,  60,  and  66  into 
eq.  58 »  one  obtains 


or  rb.  as 

66 1 

67. 


NPa  Aa 


exp(-4rrFaoarac/3kT) 


(2ifmakT)8  PaB(Pj[+PJ) 


00 

[  •xP<-WaFaB<ra-rae>a/l'T)dra 


68. 


Cancelling  terms  in  eq.  68  and  evaluating  the  integral  (  a 
standard  error  function) ,  it  is  seen  that 


2cr 


HA. 


Pa 


LtmaFaBy 


exp 


ldrro-^P, 

a  a 


p.. 


3kT(P"-P oT 


69 


where  p«=p"+p"  The  second  term  in  eq.  5 3  is  similarly  evaluated  by 
a  o 

substituting  equations  52b,  63,  and  67,  evaluating  the  integral, 
and  eliminating  rbc  by  using  eq.  64.  The  result  is 


NAb 


2c, 


* 


tfmbV' 


IL 


exp  - 


l6TTObFb 


70. 


P” 


3kT(P"-Po)‘ 


Combining  equations  69  and  70  and  expressing  o^F^  in  terms  of 
by  using  eq.  50b,  the  nucleation  rate  at  the  interface 

£1  21 

between  two  immiscible  liquids  is, 


N 


* 


* roaPaB/ 


.. — 3. - 

mf  P" 


i|Vb 

■b p" . 


exp  - 


16itc3f 
_ a  a 

3kT(P"-Po)' 


71 


where  N  =  N2/3  +  N2^3 
a  b 

The  discussion  on  pages  8  to  10  emphasized  that  the  smallest 
transfer  area  for  an  a  molecule  on  a  bubble  formed  entirely  in 
liquid  b  may  be  taken  as  the  surface  area  of  a  single  a  mole¬ 
cule.  This  suggests  that  the  transfer  areas  (i.e.,  bubble  surface 
areas)  given  by  equation  52  can  be  alternatively  expressed 
as  follows  1 

Sa  =  'WraAa.+  A.U-V  72a  > 

and 

sb  =  2fTTrbAb  +  A^1-^) 


72b. 


It  is  interesting  to  compare  the  formula  for  the  nucleation 


rate  at  a  liquid-liquid  interface  derived  in  reference  4  (re¬ 
peated  as  eq.  7^  below)  with  equations  71  and  73« 


f°aV, 


’  l6*gaFa  \ 
3kT(P"-Po)2  j 


74. 


This  equation  was  derived  assuming  that  bubble  growth  was  due 
solely  to  a  molecules  being  absorbed  on  transfer  area  S 
for  bubble  formation  as  shown  in  figs.  3c  and  6a.  The  important 
point  to  notice  in  eq,  7^  is  that  if  M  =1»  then  J=0.  But 

M  *=1  corresponds,  as  will  be  shown,  to  a  spherical  bubble 

a 

formed  in  liquid  b  and  tangent  to  the  interface.  The  fact  that 
J=0  then  implies  that  no  bubble  will  ever  form  in  this  way  (at 
least  by  the  mechanism  of  random  molecular  collisions).  Yet 
this  geometry  should  be  permissible.  This .result  .would  also 
occur  in  equations  71  and  73  by  setting  A^O  and  Ab  =0  which 
is  equivalent  to  assuming  that  no  b  molecules  are  in  a  nucleus 
formed  at  the  interface  (i.e,,  zero  transfer  area  for  b  molecules) 
A  slightly  inconsistent  recommendation  is  made  that  the 
bubble  gas  pressure,  P",  is  the  sum  of  the  partial  pressures 
of  a  and  b.  The  question,  of  course,  is  how  can  b  molecules 
enter  the  nucleus  if  its  transfer  area  is  zero?  This  is  best 
explained  by  noting  that  the  nucleation  rate  given  by  eq.  7^ 
was  derived  from  an  equation  similar  to  eq.  li 


To  produce  eq,  74  from  equations  66,  64,  52a  and  the  above 
expression,  p  must  be  given  by  the  following  (if  P"=P”+P" ) t 

21  O 

i  fi  -  - ^ - “X 

(2frm  kT)? 

3. 


which  does  not  seem  quite  right. 

The  problem  is  that  there  are  two  transfer  areas  instead 
of  one  and.  that  eq.'l  should  have  been  used  to  determine  the 
nucleation  rate.  One  can  argue,  however,  that  these  "errors" 
have  a  small  effect  in  that  they  appear  only  in  the  pre -exponen¬ 
tial  term  in  the  expression  for  the  nucleation  rate.  Neverthe¬ 
less  it  will  be  shown  that  the  absurd  result  J=Q  can  occur 
under  certain  limiting  conditions  in  eq.  74. 

The  particular  conditions  which  render  equation  71  or 
equation  73  valid  were  determined  in  reference  4  (using  eq.  74) 
in  a  rather  elegant  manner.  Combing  equations  46a,  47,  and  55 
together  with  the  identity  sin(e)  =  (l-cos(e))2,  the  character¬ 
istic  lens  dimension,  a  in  fig,  6a, can  be  expressed  as  followsi 


This  formula  was  factored  in  reference  4  and  appears  as  eq.  75 


Prom  eq.  75  "the  conditions  for  determining  the  location  of  bub¬ 
ble  formation  at  a  liquid-liquid  interface  shown  in  fig.  3 


are  now  seen  to  depend  only  on  the  appropriate  surface  and 
interfacial  tensions.  Depending  on  these  values  the  following 
cases  may  mathematically  occur « 


1)  a  is  real i 

2)  a  is  zero}  and 

3)  a  is  imaginary. 

Depending  on  whether  the  second  term  in  brackets  in  eq.  74  is 
positive,  zero,  or  negative  gives  the  required  conditions 
because  physically  <rat><cr^L+ab  always  occurs, 
a  is  real  only  when 


°ab  >  IV'bl 

The  lens  radius  is  zero  when 


77 

78 
79. 


80 

81. 


These  three  conditions  are  discussed  below. 

1)  a  is  real> 

In  this  case  equation  76  applies.  A  bubble  exists  between 
liquids  a  and  b  as  illustrated  in  fig.  3o,  and  equation  71  or 


equation  73  gives  the  expression  for  the  nucleation  rate, 
2)  a  is  zero i 

Either  of  equations  78  or  79  apply  in  this  case. 

i .  &  -  ar  +cr  71 


When  equation  78  applies  it  can  easily  be  shown  that 
(see  equations  46,  50,  and  52) 

P  85 

a  cr  a 

Fb=l  85 

A  =0  .85 

a 


V* 

Vs-1 


82d , 


82f. 


Under  these  conditions  equation  73  reduces  to  the  following 


equation* 


b\rb 

/P" 


(128rrB<^): 


exp  - 


l6^°b  - 

3kT(P"-Po) 

eq.  83 


Equation  71  with  A*=0  becomes 

/o.  \  i  DM 


J  =  N  $v) 


p-  / 

lb  exp  -  - , 

P"  \3kT(P”-Po)‘ 


Substituting  equation  82  into  eq.  74  gives,  in  addition,  J^O.. 
Comparing  equations  71  and  73  with  74  shows  why  this  occurs. 
There  is  no  transfer  area  for  b  molecules  so  that  there  is  only 
one  pre -exponential  term  in  eq,  74  and  this  term  is  zero  when 


equation  78  is  satisfied. 

It  may  be  recognized  that  equation  83b is  the  formula  for 
homogeneous  nucleation  in  liquid  b  for  bubble  growth  as  illustrated 
in  fig.  3d  (except  N  is  now  the  number  density  of  surface  mole¬ 
cules).  This  is  not  surprising  because  Fb=l  and  hence 
V-=  4/3 IT r^  so  lhat  the  nucleus  is  spherical  of  radius  r^. 

The  work  to  form  the  critical  size  nucleus  is,  from  the 
exponential  term  in  equation  83o-(or  83l) , 


84. 


Equation  84  is  exactly  the  formula  which  would  be  obtained  for 
the  work  to  form  a  spherical  bubble  of  critical  size  in  the  bulk 
of  liquid  b.  At  the  condition  under  present  discussion,  the 
bubble  forms  at  the  interface  and  is  tangent  to  the  surface 
because  9  =  0  and  (J>  =180  when  Ma=l  and  Mb=-1 . 


11  •  %  "  °ab+ffa 

79 

When  eq.  79  is  satisfied,  it  can  be  shown  that 

P  =1 
a 

85a 

3* 

it 

pj 

& 

85b 

Aa=1 

85c 

0 

11 

85d 

H 

II 

«« 

85e 

rH 

II 

85f 

Under  the  conditions  given  by  equations  85»  eq.  73  reduces  to, 

)  p,  +  Mb( f-vVs 


8? 


When  At,  =0,  equation  71  reduces  to 


16  ir<f? 
_ a 


3kT(P"-Por 


The  second  term  in. brackets  in  eq.  86  is  expected  to  be  negli¬ 
gible  compared  with  the  first  term  because  P£  »  P£.  Eq,.  86 
then  further  reduces  to  eq.  87. 

The  work  to  form  the  critical  size  bubble  is 


W  = 


16ffo; 


3.(PM-Por 


88. 


Equations  87  and  88  are  exactly  what  one  would  obtain  for 
homogeneous  nucleation  in  the  bulk  of  liquid  a  except  that  N 
would  be  the  number  density  of  a  molecules  and  not  the  number 
density  of  molecules  at  the  interface  between  the  two  liquids. 
In  addition  for  homogeneous  nucleation  in  the  bulk  of  liquid  a, 
only  a  molecules  would  be  in  the  bubble  so  that  P“/P"  =  1  and 

""  a 

P£/P"  =  0. 


A  spherical  bubble  is  formed  when  eq.  79  is  satisfied  because 
F  =1  in  eq.  50b.  In  addition  0  =0  and  9=180  (from  eqs.  46a  and 

fl. 

46b)  so  that  the  bubble  is  tangent  to  the  interface  and  sur¬ 
rounded  by  liquid  a  as  illustrated  in  fig.  3b. 

3)  a  is  imaginary! 


This  is  the  most  interesting  situation  because  there  is  no 
rigorous  mathematical  limit  to  which  eq.  71  or  eq.  73  approaches. 
Moreover  this  condition  is  requently  encountered  in  nature. 

Simply  stated,  when  a  is  imaginary,  a  nucleus  cannot  form 
anywhere  at  the  interface  and  must  grow  away  from  the  inter- 


face  in  the  bulk  of  either  liquids  a  or  b  (figs.  3a  and  3e). 
Equations  80  and  81  are  discussed  separately  below. 

ga  *  qb'l~°ab  80> 

When  equation  80  is  satisfied,  it  can  be  showft  that 

Ma  >  1. 

“b  <-1  • 

and  >  1# 

But  the  cosine  of  an  angle  can  never  fee  greater  than  1.  Therefore 
there  is  no  bubble  configuration  at  the  inteiface  describable 
by  the  angles  0  and  ,  and  hence  no  bubble  exists  at  the  inter¬ 
face  (at  least  in  a  mathematical  sense). 

The  fact  that  >  1  implies  that  if  a  bubble  could  grow 
at  the  interface,  the  energy  of  its  formation  to  the  critical 
size  would  be  less  in  the  bulk  of  liquid  b  than  at  the  interface. 


This  is  so  because  from  equations  55  and  67  at  r=rtc 

w  „ 

3(P"-P0)2 


89, 


while  the  work  to  form  a  bubble  of  critical  size  in  the  bulk  of 


liquid  b  is  given  by  eq,  84  which  is  smaller  than  eq.  89  by 

r.  In  addition  the  nucleation  rate  is  greater  by  a  factor  of 

b 

10^  when  a  nucleus  forms  homogeneously  within  the  bulk  of  a 
liquid.  Both  of  these  factors  tend  to  increase  the  probability 
that  nucleation  will  occur  in  the  bulk  of  liquid  b  rather  that 
at  the  interface.  The  choice  is  essentially  between  the  loca¬ 


tions  illustrated  in  figs.  3c,  3d,  and  3e. 

For  homogeneous  nucleation  in  the  bulk  of  liquid  b  there  are 


no  a  molecules  in  the  bubble  so  that  P"=0  and  ^  =0.  The  bubble  is 


spherical  which  implies  that  A-=l.  Equation  73  then  becomes 


N*te) 


£  /  16  fra? 


exp  - 


3kT(P^-Po)‘ 


90. 


Equation  80  is  most  likely  to  be  satisfied  by  water-fuel 

emulsions  stabilized  by  an  emulsifying  agent.  For  such  emulsions 

P"  »  P"  (liquid  a  is  water  and  b  is  fuel).  Therefore,  the  ar- 
a  o  ~  *” 

guments  presented  on  pages  8  to  10  apply  in  this  case  as  measured 
superheat  temperatures  are  just  too  low  to  be  explained  by 
homogeneous  nucleation  in  liquid  b.  This  means  there  must  be  a 
non-zero  transfer  area  for  a  molecules.  Assuming  that  the  bubble 
forms  close  enough  to  the  interface  so  that  the  dominant  growth 
mechanism  is  that  of  random  molecular  absorption  on.  the  bubble  surface 
(i.e,,  diffusion  is  neglected),  the  nucleation  rate  may  be 
approximated  by  equation  83a.  The  suggestions  are  then  the 
following*  When  P"  »  P£,  equation  83b  or  equation  83b(AcfO  in 
eq.83a)  will  apply.  If,  however,  P"  «  P",  the  expected  result 

CL  D 

is  equation  90  (P"=0  and  Av=l  in  eq.  73).  Experimental  results 

a  D 

are  in  quite  good  agreement  with  eq,  83a  (P"  »  P”)  lending  sup- 

a  D 

port  t  -.  these  recommendations, 

ii.  <r.  >  <T  +<T  .  81 

b _ a  ab 

When  this  condition  is  satisfied,  it  can  be  shown  that 

M  <-l 
a 

v1 

v1- 

When  F_>1  it  takes  less  work  to  form  a  spherical  bubble  of  critical 


size  in  the  hulk  of  liquid  a  than  at  the  interface.  This  is  so 
because  for  a  bubble  dividing  the  interface  (fig.  3c) 

w  ,  91. 

'  3(P"-P0)2 

which  is  larger  than  eq.  88  by  F  ,  The  use  of  eq.  88  sug- 

a 

gests  that  F  =1  (spherical  bubble),  Aa=l,  and  P£=0  (no  b  mole¬ 
cules  in  the  nucleus)  for  a  bubble  forming  away  from  the  inter 
face  in  the  bulk  of  liquid  a.  Substituting  these  conditions 


into  eq.  73,  it  is  seen  that 

/  \  ^ 


J  ■ 


exp  - 


1 6?rof 

_  a 

3kT(P"-P  ) 
a  o 


If  p«  «  p",  the  same  dilemma  will  arise  as  discussed 
a  b 

before i  the  bubble  must  grow  close  enough  to  the  interface  for 

b  molecules  to  enter.  This  means  that  the  nucleus  must  be  in 

contact  with  the  interface.  The  above  arguments  can  be  used  to 

show  that  equation  86  will  approximate  the  nucleation  rate  in 

this  case.  This  suggests  that  when  P£  »  P"  and  o.  >  o'  +o  . 

b  a  b  —  a  ab, 

equation  86  can  be  used  to  approximate  the  nucleation  rate. 


4.  Summary  of  Theoretical  Discussion 


Kfl«ll 


The  purpose  of  the  experimental  program  was  to  determine 
which  of  equations  76  to  81  was  satisfied  for  the  particular 
water-fuel  emulsions  used  here,  and  whether  the  corresponding 
expression  for  the  nucleation  rate  could  be  used  to  predict  the 
temperature  at  which  nucleation  would  occur  during  burning  of 
the  emulsified  fuel  droplet. 

The  ‘  method  used  here  consisted  of  preparing  an  emulsion 
of  two  liquids,  a  and  b,  and  then  heating  a  sample  drop.  Such 
a  drop  is  shown  schematically  in  fig.  7a.  The  interface  of  the 
internal  phase  and  surrounding  liquid  can  be  viewed  as  a  flat 
interface  (illustrated  in  fig.  7b).  It  is  assumed  that  the  tiny 
droplets  of  liquid  a  are  far  enough  apart  so  that  the  surrounding 
liquid  can  be  considered  as  an  infinite  mass.  The  geometry  of 
fig.  7b  then  becomes  that  of  fig.  3. 

For  an  emulsion  drop  heated  in  air,  the  maximum  temperature 
will  be  limited  by  the  boiling  point  of  liquid  b.  Unless  the 
nucleation  temperature  at  the  interface  between  the  two  liquids 
is  lower  than  this  temperature,  nucleation  will  not  occur.  The 
b  liquid  may  also  have  to  be  superheated.  This  suggests  the 
classic  nucleation  experiment  in  which  a  drop  of  a  test  liquid 
is  injected  in  the  bottom  of  a  column  filled  with  a  heavier 
immiscible  liquid  under  a  temperature  gradient.  The  rising  drop 
is  progressively  heated  until  it  explodes.  The  maximum  droplet 
temperature  attainable  in  this  experiment  is  the  limit  of  super- 


heat  of  liquid  b  (assuming  the  boiling  point  of  the  column  fluid 
is  higher  than  this  temperature).  Since  the  limit  of  superheat 
of  a  liquid  is  usually  around  90 %  of  its  critical  temper¬ 
ature!  considerable  superheat  of  the  liquid  can  therefore  be 
sustained  in  the  emulsion.  One  is  not  limited  by  the  boiling 
point  of  liquid  b  in  this  experiment. 

The  design  of  the  apparatus  is  basically  that  described  in 
reference  11. 

2.  Experimental  Apparatus  and  Method 

A  schematic  illustration  of  the  apparatus  is  shown  in  fig.  8. 
The  column  consisted  of  a  pyrex  tube  4.2cm  inside  diameter  and 
68cm  long.  The  central  37cm  was  wound  with  32  turns  of  nichrome 
wire  and  cemented  in  place  with  Sauereisen  cement.  The  spacing 
between  turns  varied  from  2cm  at  the  lower  end  to  ,5cm  at  the 
upper  end.  Power  to  the  nichrome  wire  was  provided  by  an  AC 
voltage  supply  controlled  by  a  variac.  Temperature  was  measured 
by  one  #28  gage  copper-constantan  thermocouple  enclosed  in  a 
“L"  shaped  pyrex  tubing  jacket  inserted  directly  into  the  col¬ 
umn  heating  liquid.  The  thermocouple  was  connected  to  a  Leeds  & 
Northrup  Model  8696  potentiometer  and  calibrated  to  about  .1C. 

A  #10  rubber  stopper  was  inserted  irv+o  eac  and  of  the  column. 

The  thermocouple  tubing  jacket  and  a  condenser  were  fitted  into 

the  top  robber  stopper.  The  condenser  reduced  the  amount  of 

column  heating  liquid  lost  by  vapor  generation  during  heating 

and  prevented  harmful  vapors  from  entering  the  laboratory.  A 
-nitrogen  gas  atmosphere  above  the  heating  liquid  was  provided 


■to  reduce  oxidation  and  subsequent  darkening  of  the  column 
heating  liquid  (glycerine  in  this  case).  Efforts  were  made  to 
insure  that  atmospheric  conditions  were  maintained  by  only  lightly 
purging  the  gas  above  the  heating  liquid. 

Droplet  injection  was  by  a  hypodermic  needle.  A  1/8  inch 

diameter  by  1.5  inch  long  brass  tube  with  a  3/4  inch  diameter 

toy  3/4  inch  long  brass  cylinder  soldered  on  the  end  was  fitted 

into  the  center  of  the  bottom  rubber  stopper.  A  rubber  septum 

was  placed  over  the  brass  cylinder  which  could  be  punctured  by 

the  needle  (see  fig.  8).  A  number  22  gage  hypodermic  needle 
2  inches  long  fitted  with  a  lcc  glass  syring  was  used  to  inject 
the  drops. 

Because  of  the  necessity  of  knowing  accurately  fuel  proper¬ 
ties,  only,  pure  fuels  were  used  in  the  emulsions.  Test  droplets 
consisted  of  water-in-fuel  emulsions.  The  emulsions  were  pre¬ 
pared  by  pouring  a  selected  volume  percentage  of  water  in  the 
fuel  and  mixing  by  a  counter  rotating  propeller  (fig.  9).  Be¬ 
cause  of  the  high  water-fuel  interfacial  tension  (typically 
between  30  dynes/cm  and  60  dynes/cm  at  20C),  it  was  necessary 
to  stabilize  the  emulsion  by  adding  an  emulsifying  agent.  Other¬ 
wise  the  water  would  coalesce  and  collect  in  the  bottom  of  the 
mixing  beaker*  The  emulsifying  agent  was  first  mixed  in  the  fuel 
and  the  water  added  later.  Various  samples  containing  between 
1%  and  2%  emulsifying  agent  were  prepared.  The  emulsifying  agent 
was  a  mixture  of  60fo  Span  85  and  40#  Tween  85  (available  from 
ICI  America,  Inc.),  This  particular  mixture  was  experimentally 
found  to  stabilize  higher  alkane-water  emulsions  for  at  least 


15  minutes  at  room  temperature  which  was  ample  time  to  draw  out 
a  freshly  made  sample  and  inject  a  drop  into  the  heating  column. 
For  other  water-fuel  combinations  this  particular  emulsifying 
agent  proved  unsatisfactory.  In  particular  a  methanol-water 
emulsion  could  not  be  prepared  in  this  way  without  the  water 
coalescing. 

The  emulsions  were  made  up  of  83#  fuel,  15#  water  and  2% 
emulsifying  agent  (  by  volume).  The  four  water-fuel  emulsions 
tested  were  the  following i 

1)  water-decane j 

2)  water-dodecane t 

3)  water-tetradecanej  and 

4)  water-hexadecane. 

A  test  run  proceeded  by  stabilizing  a  temperature  gradient 
in  the  column  and  injecting  a  drop  of  the  test  emulsion.  Meaure- 
ments  consisted  of  recording  the  height  of  explosive  boiling  and 
the  corresponding  temperature.  Explosions  of  large  drops  (>  2mm 
diameter)  occasionally  upset  the  column  temperature  gradient, 
and  about  15  minutes  to  30  minutes  was  required  to  re-stabilize 
the  gradient.  Fig.  10  shows  a  typical  temperature  gradient  in  the 
column  (the  temperature  was  adjusted  to  286c  at  the  top  for 

a  pure  decane  experiment).  The  temperature  gradient  was  very 
flat  at  the  top  of  the  column,  varying  about  2C  per  10cm. 

Whenever  possible,  the  gradient  was  adjusted  to  cause  explosions 
in  thia  region. 

No  liquid  is  currently  available  which  is  immiscible  both  with 
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water  and  the  fuels  listed  on  the  previous  page,  heavier  than 
both,  and  with  a  boiling  point  greater  then  about  300C.  The 
best  liquid  for  use  as  a  heating  medium  appeared  to  be  glycerine 
which,  unfortunately,  is  very  soluble  in  water.  A  certain 

i 

amount  of  surface  dissolution  of  the  water  into  the  glycerine 
therefore  had  to  be  allowed  to  obtain  nucleation  temperatures 
by  the  method  used  here.  Nevertheless,  if  it  is  assumed  that 
nucleation  occurs  at  a  water-fuel  int  .rface  within  the  drop, 
the  corresponding  temperature  should  be  predictable  from  ’  e 
theory  presented  here.  The  assumpt  on  is  being  made  tv  a 
pure  fuel  shell  surrounds  the  intei  lal  emulsion  and  no  water 

can  penetrate  the  fuel  glycerine  i  .terface. 

Additional  experiments  were  i  un  using  pure  hydrocarbon  drops 
to  test  existing  theories  and  to  oecome  familiar  with  the  op¬ 
eration  of  the  apparatus.  The  hydrocarbons  were  obtained 

from  Humphrey  Chemical  Co.  and  cf  99.6%  purity  or  better.  The 
following  pure  liquids  were  useii  n-pentanej  n -hexane j  benzene ; 
n-octanej  and  n-decane. 


3.  Experimental  Results 

Tables  I  and  II  list  the  experimentally  measured  superheat 
temperatures  for  the  pure  hu  irocarbons  and  emulsions  tested. 

The  nucleation  temperatures  of  the  pure  hydrocarbons  are 
completely  consistent  with  'alues  reported  in  the  literature.  In 
all  cases  the  injected  droplets  exploded  with  a  sharp  report, 
the  largest  droplets  (<^-,2mf  diameter)  actually  shaking  the 
apparatus . 


Droplet  explosions  for  the  pure  hydrocarbons  occured  in  a 
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very  narrow  temperature  range.  Stating  an  actual  mumerical  bound 
is  somewhat  fortuitous.  As  an  example,  for  the  pure  n-pentane 
drops  the  thermocouple  was  positioned  at  the  observed  level  of 
explosion.  ,For  all  injected  n-pentane  drops,  this  level  did 
not  have  to  be  changed  because  all  drops  exploded  at  the  same 
height. 

The  temperatures  listed  in  Table  I  represent  an  average  of 
the  recorded  values  and  are  accurate  to  within  at  least  ±  1C. 

The  superheat  temperature  of  n-decane  was  somewhat  diffi¬ 
cult  to  obtain.  The  measured  value  of  about  283C  is  only  7C 
lower  than  the  boiling  point  of  glycerine.  At  this  high 
temperature,  stress  cracks  in  the  glass  column  were  observed,  and 
on  one  occasion  the  bottom  of  the  pyrex  tube  literally  broke  off. 
The  glycerine  could  not  be  maintained  for  more  than  £  hour 
at  283C  before  such  trouble  began  to  oocur.  However,  once  the 
proper  gradient  was  established,  and  working  swiftly,  the  super¬ 
heat  temperature  of  n-decane  exhibited  the  same  narrow  range 
of  explosive  vaporization  as  the  other  hydrocarbons  tested. 

Higher  alkanes  such  as  dodecane,  hexadecane,  and  tetradecane  have 
theoretical  limits  of  superheat  well  in  excess  of  300C  which 
just  cannot  be  measured  using  glycerine  as  the  heating  medium. 

One  would  also  be  better  off  using  a  quartz  column  to  support 
such  temperatures.  The  highest  superheat  temperature  ever 
measured  at  latm  was  287C  for  cyclo-octane  (reference  11).  Need¬ 
less  to  say  the  apparatus  must  have  pushed  to  its  limit  to  measure 
this  temperature. 

Injected  droplet  size  was  determined  by  measuring  the  rate 


rise  of  a  selected  drop  over  a  specified  distance  at  a  given 
temperature.  The  low  Reynolds  number  solution  of  the  terminal 
velocity  of  a  liquid  drop  rising  freely  in  an  immiscible  heavier 
liquid  as  derived  in  reference  12  was  then  used  to  correlate  the 
measured  rate  of  rise  with  the  drop  diameter.  Droplet  diameters 
for  both  the  pure  fuels  and  emulsions  ranged  between  .  5®m 
(500  microns)  and  about  2mm  (2000  microns).  For  droplets  larger 
than  about  2mm  diameter  thermal  equilibrium  cannot  be  assumed 
(ref.  13).  Data  from  drops  of  this  size  are  not  included  in 
Tables  I  and  II. 

Internal  dispersed  phase  (water)  size  of  the  emulsions  was 
difficult  to  determine.  Examination  of  a  typical  emulsion  sample 
under  a  microscope  at  room  temperature  revealed  some  coalescence 
and  internal  motion.  Generally,  though,  the  order  of  magnitude 
of  *he  internal  water  drop' diameter*  appeared  to  range  from  about 
.01mm  1  to. .1mm. .  For  even  the  smallest  drops  the  assumption  of 
a  flat  plane  interface  for  nucleation  is  a  good  one  because 
the  diameter  ofa spherical  nucle  of  critical  size  is  about 
10“^cm,  and  10“^  »  10~^.  Techniques  for  varying  the  internal 
phase  size  of  emulsions,  as  well  as  photomicrographs  of 
typical  water-fuel  emulsions,  are  given  in  reference  14. 

Rising  emulsion,  drops  left  a  visible  trail  behind  them. 

This  indicated  that  water  was  slowly  dissolving  in  the  glycerine. 
Some  settling  of  the  water  was  also  observed.  As  a  droplet  rose, 
the  formation  of  a  visible  interface  between  the  clear  fuel  and 
the  milky  internal  emulsion  was  seen.  An  illustration  of  this  is 


shown  in  fig.  11.  Whether  this  interface  was  formed  as  a  result 
of  drainage  or  settling  of  the  heavier  water,  or  a  combination 
of  both,  is  not  known.  In  any  case  this  was  more  pronounced  in 
the  smaller  droplets  (i.e,,  less  than  around  ,5mm  diameter). 

For  this  size  drop#water  dissolution  was  so  pronounced  that  be¬ 
fore  the  droplet  even  reached  its  nucleation  temperature  it 
had  become  completely  clear  indicating  that  all  the  water  had 
dissolved  at  the  surface.  The  resulting  pure  fuel  drop  would 
then  rise  right  to  the  top  of  the  column  without  exploding. 

It  was,  therefore,  necessary  to  have  the  injected  droplet  large 
enough  so  that  there  would  still  be  some  of  the  internal  emulsion 
left  by  the  time  the  droplet  reached  its  nucleation  temperature. 
In  this  regard,  the  size  of  the  injected  droplet  rather  than 
the  percentage  of  water  in  the  emulsion  is  more  important. 

Emulsions  with  7 %  water  and  1 %  emulsifying  agent  were  also 
tested.  No  change  in  the  maximum  observed  nucleation  temperature 
was  detected.  There  also  appeared  to  be  no  dependence  of  the 
nucleation  temperature  on  injected  droplet  diameter.  The  wide 
variety  of  internal  water  sizes  described  above  also  provide 
evidence  of  the  independence  of  the  nucleation  temperature  on 
.droplet  ..r=  size.  These  observations  are  consistent  with  the 
theory  presented  before  and  the  data  reported  in  the  literature. 

Emulsion  explosions  were  much  less  violent  than  pure  fuel 
explosions.  Rather  than  a  sharp  report,  an  almost  inaudible 
"puff”  accompanied  each  explosion. 

Fully  90%  of  the  injected  emulsion  droplets  were  observed  to 
have  air  bubbles  attached  to  them.  This  was  due  to  the  counter 


rotating  propeller  causing  significant  cavitation  in  the  fuel 
in  the  mixing  beaker.  Expansion  of  the  air  bubbles  on  the  rising 
drops  could  be  visually  followed. These  expansions  occured  within 
the  first  5cn>  of  the  heated  section  of  the  column  and  were  thus 
easy  to  detect.  The  column  temperature  gradient  was  temporarily 
upset  by  this  ani usually  required  15  minutes  to  re-stabilize. 

Even  without  attached  air  bubbles,  droplet  explosions  ..  .of 
the  fuel  emulsions  appeared  to  occur  over  a  much  wider  temperature 
range  than  the  pure  .  liquids.  Imperfectly  wetted  foreign 
particles  in  particular  can  lower  superheat  temperatures.  Table  II 
lists  the  highest  recorded  temperature  for  the  particular; 
water-fuel  emulsion  tested  (in  recording  only  the  highest  measured 
temperature,  the  reader  is  reminded  of  the  work  described  in 
references  13  and  15). 

IV.  Discussion  of  Results 

1.  Surface  Tension 

By  determining  conditions  which  would  result  in  real  and 
imaginary  solutions  to  the  equation  for  the  "radius"  of  a  vapor 
lens  at  a  liquid-liquid  interface  -  eq.  75  -  criteria  were  for¬ 
mulated  for  determining  where  a  nucleus  of  critical  size  would  ap¬ 
pear.  These  criteria  were  found  to  depend  only  on  the  surface 
and  interfacial  tensions  of  the  mutually  saturated  liquids. 

Surface  tension  data  of  mutually  saturated  liquids  generally 


do  not  exist.  Most  data  reported  in  the  literature  are  for  liquids 
saturated  either  with  their  own  vapor  ,  or  in  air.  The  assumption 
will  be  made  that  the  error  in  using  the  pure  liquid  component 


surface  tensions,  o'  and  o',  ,  is  small. 

a  o 

The  surface  and  interfacial  tension  ■sariation  with  temper¬ 
ature  must  be  known  to  use  any  of  the  for  aulas  for  the  nucleation 
rate  previously  derived.  Because  these  c ata  are  usually  not 

I 

available  at  liquid  superheat  temperatur  .*s,  it  is  necessary  to 
extrapolate  to  these  temperatures.  The  error  incurred  by  this 
extrapolation  accounts  for  the  dominant  error  in  the  prediction 
of  nucleation  temperatures. 

Pure  liquid  surface  tension  valuer  covering  a  wide  variety  of 
liquids  appear  in  reference  16.  Data  for  all  liquids  reported 
are  given  in  the  range  of  about  20C  tj  150C.  For  extrapolation 
to  higher  temperatures,  it  is  useful  to  correlate  these  data 
with  the  following  formula  1 

cr  =  cr0(l-T/Tc)n  93 

where  T  is  the  critical  temperatu  ’e  and  a  and  n  are  constants 
c  o 

whoso  numerical  values  depend  on  the  particular  liquid.  Correlations 

in  the  form  of  eq.  93  for  _  all  liquids  studied  here,  except 

benzene  and  water,  were  found  in  the  literature.  For  benzene 

and  water  a  computer  library  sub  ’outine  was  used  to  correlate 

the  data  given  in  reference  17  (for  water)  to  a  least  squares 

curve.  Values  of  the  critical  ;emperature  were  taken  from 

reference  18.  Appendix  I  lists:  the  constants  in  eq.  93  and 
the  least  squares  constants  of  the  equation 

^  s  g  ip  +  a  ^  ...  9^ 

o  . .  2 

for  benzene  and  water. 


Interfacial  tension  data  at  high  temperature  were  virtually 
non  existant  except  for  decane-water.  The  interfacial  tension 

•for 

at  only  20C  could  be  located/^  For  hexadecane-water  and  dodecane- 
water  no  values  at  any  temperature  could  be  found.  The  extent 
of  the  interfacial  tension  iata  for  the  emulsions  tested  is 
summarized  in  Table  III  below. 

Table  III 

Interfacial  Tension 

Water-Tetradecane i  crab  =  52.2  dynes/cm  (@20C)  (ref.  19) 
Water-Decane i  0&b=  50.066+.0027247Po-.1205(T-298.15)  (ref. 20) 

A  formula  relating  the  inter facial  tension  to  the  pure 
liquid  component  surface  tensions  was  derived  in  ref.  21* 

°ab  *  °a  +  ab  -  2<<40b^  95- 

The  last  term  in  eq.  95  represents  the  dispersion  force  com¬ 
ponent  of  the  interfacial  tension.  For  the  hydrocarbons  used 
here  this  is  equal  to  the  pure  liquid  surface  tension  (ref.  21) 
so  that  0^=0^.  The  dispersion  component  of  water,  <?a,  was 
estimated  in  reference  9  to  be 

where  K  is  a  constant  and  J and  are  liquid  and  vapor  den¬ 
sities  of  water  respectively.  Combining  this  with  eq.  95  gives 
the  relation  between  surface  and  interfacial  tension  at  a 
water-fuel  interface  and  its  subsequent  variation  with  temperature i 


where  the  subscript  "o"  denotes  the  reference  temperature  at 

which  <j  ,  is  known, 
ao 

For  a  tetradecane -water  interface  at  25c,  <J  =71.965dyne/cm, 

£L 

tfb=26.126dyne/cm,  and  crat)=52.2dyne/cm.  From  eq.  97  (and  thermo¬ 
dynamic  tables  fbr  water  properties),  it  is  found  that 
12  3  2 

K  =  20.401  cm  /g^-sec  .  Similarly  for  a  water-decane  interface 
it  can  be  shown  that  K  =  21.79  cm  /g-'-sec  ,  and  for  octane- 
water,  K  =  22.8  cml2/g-^-sec2 ,  These  results  suggest  that  the 
interfacial  tensions  of  alkane-water  systems  do  not  differ 
greatly. 

The  location  of  bubble  formation  can  now  be  determined  from 
the  conditions  summarized  on  pages  43  to  45,  Consider  as  an 
example  a  water-tetradecane  emulsion.  The  highest  measured  nucleation 
temperature  is  about  260C  (Table  II).  At  this  temperature 
or  *  23.85  dyne/cm,  o'.  «7.9699dyne/cm,  and  from  eq.  96  with 

cl  D 

Ke20.401cm12/g^-sec2,  cat)s*17.078dyne/cm.  Since 

17,078  >  |  23.85-7.9699|=  15.88, 

eq,  76  is  satisfied.  This  means  that  one  should  expect  a  bubble 
to  form  between  a  water-tetradecane  interface  as  illustrated 
in  fig,  3c.  Equation  73  or  eq.  71  then  provides  the  expression 


for  the  nucleation  rate.  Similarly  for  a  decane-water  interface 
at  230C,  tf  «30 . 97dyne/ cm ,  <Jb«5.787dyne/cm,  and  from  Table  III 
oab«>25.366dyne/cm.  Since 

'25.366  >  |30.97  -  5.78?|=  25.183, 

eq.  76  is  again  satisfied  (though  just  barely). 

It  is  interesting  to  note  that  equations  76  to  81  are  also  the 
conditions  for  determining  whether  one  liquid  will  spread  on 
another.  A  detailed  discussion  of  this  will  not  be  given  here 
(see  ref.  22).  Briefly  when  eq,  76  is  satisfied,  a  drop  of 
liquid  b  on  liquid  a  will  form  a  lens.  Equations  80  and  81,  of 
which  equations  78  and  79  are  limiting  cases,  are  the  conditions 
of  spreading  of  one  liquid  on  another.  For  example,  when  eq.  80 
is  satisfied,  a  drop  of  liquid  b  on  liquid  a  will  immediately 
spread  on  the  surface  of  a.  This  is  equivalent  to  a  drop  of  b 
resting  on  the  surface  of  a  with  a  zero  contact  angle,  A  similar 
explanation  applies  to  eq.  81  for  a  drop  of  liquid  a  on  liquid 
b.  The  implication  is  that  when  a  surface  is  completely  wetted 
by  another  liquid  so  that  equations  80  or  81  applies,  a  spherical 
bubble  will  form  either  in  the  bulk  of  the  wetting  liquid  (fig,3e) 
or  at  the  interface  surrounded  by  the  wetting  liquid  (fig.  3d), 
depending  on  which  liquid  is  more  volatile. 

For  a  water-decane  interface  at  20C,  the  data  in  Table  III 
and  Appendix  I  indicate  that  eq.  76  is  satisfied.  From  the 
above  discussion  this  means  Hi  at  a  drop  of  decane  will  rest 
as  a  lens  on  the  water,  which  is  easily  verified.  This  can  also 
be  verified  for  the  other  water-fuel  systems  in  Table  II, 


When  2%  by  volume  of  emulsifying  agent  is  mixed  with  any 
of  the  pure  hydrocarbons  listed  in  Table  II,  and  a  drop  of  this 
hydrocarbon-emulsifying  agent  mixture  placed  on  water,  the  drop 
spreads  instantly  on  the  water  surface.  Eq.  76  is  no  longer 
satisfied  and  the  spreading  condition,  eq.  80,  describes  the 
water-fuel  interface.  The  stabilizing  effect  surfactants  have 
on  emulsions  is  due  primarily  to  a  lowering  of  the  interfacial 
tension  (the  emulsifying  agent  being  absorbed  on  the  water-fuel 
interface).  Although  no  measurements  were  made  here  of  the 
interfacial  tension  at  a  water- fuel  interface,  it  is  apparent 
that  if  such  data  were  taken,  very  low  interfacial  tensions 
would  be  measured. 

Depending  on  the  surfactant,  it  is  possible  to  lower  the 
interfacial  tension  from  a  typical  value  of  50dyne/cm  to 
less  than  ldyne/cm  at  20C,  Such  a  dramatic  decrease  of  inter- 
facial  tension  will  change  the  location  of  bubble  growth  from 
a  lens  forming  between  the  water-fuel  interface  (fig.  3c)  to 
spherical  bubble  formation  at  the  interface  or  in  the  bulk  of  the 
fuel.  When  this  occurs,  it  is  not  necessary  to  actually  know  a 
numerical  value  of  the  interfacial  tension.  This  is  so  because 
the  work  to  form  a  spherical  bubble  depends  only  on  the  pure 
component  surface  tension  of  the  liquid  in  which  the  bubble  grows 
(i.e.,  equations  84,  88,  and  89), 

To  determine  the  nucleation  rate  of  an  emulsion  stabilized 
by  a  surfactant,  it  is  only  necessary  to  know  which  liquid  spreads 
on  the  other,  a  condition  which  is  sure  to  be  satisfied.  The 


the  interfacial  tension  causes  a  change  in  the  nucleation  mechanism 
to  the  extent  that  the  corresponding  equation  for  the  nucleation 
rate  is  independent  of  the  interfacial  tension. 

Surfactant  absorption  at  a  water-fuel  interface  might  sug¬ 
gest  that  the  pure  liquid  component  surface  tension  does  not 
change  appreciably.  This  hypothesis  was  tested  by  measuring 
the  limit  of  superheat  of  n-pentane  and  n-decane  each  mixed  with 
2j6  emulsifying  agent  (60%  Span  85  and  40$  Tween  85).  Superheat 
temperatures  were  nearly  the  same  as  for  the  pure  liquids,  ranging 
between  145C-147C  for  n-pentane  and  280C-285C  for  n-decane. 

Given  a  nucleation  rate  and  temperature,  eq.  90  can  be  solved 
for  surface  tension.  Because  the  measured  superheat  temperatures 
were  nearly  the  same  as  that  of  the  pure  fuels,  the  pure  fuel 

+e«s*o* 

surfaceAis  essentially  unaffected  by  the  addition  of  the  sur¬ 
factant  used  in  this  study. 

The  choice  of  a  heating  medium  depends  on  whether  eq.  81 

is  satisfied.  The  hydrocarbons  used  here  spread  on  glycerine 

indicating  that  eq.  81  is  indeed  satisfied.  In  addition  since 

P"  »  P" ,  nucleation  will  occur  within  the  bulk  of  the  injected 
b  g 

fuel  drop  (homogeneous  nucleation),  and  not  at  the  fuel  glycerine 
interface.  The  results  in  Table  I  verify  this.  A  rising  emulsion 
drop  can  therefore  be  modelled  as  a  droplet  burning  in  air  ex¬ 
cept  that  now  the  droplet  temperature  is  able  to  reach  its  limit 
of  superheat  and  not  just  its  boiling  point.  The  advantages  of 
this  are  obvious  particularly  in  light  of  the  results  shown  in 
Table  II. 


2,  Vapor  Pressure 

The  component  gas  partial  pressure  in  a  nucleus  of  critical 
size  is  not  its  equilibrium  vapor  pressure,  but  is  slightly  less. 
This  is  due  both  to  bubble  curvature  and  the  fact  that  the  ex¬ 
ternal  liquid  pressure  is  at  the  ambient  rather  than  the  com¬ 
ponent  equilibrium  vapor  pressure. 

A  nucleus  of  critical  size  is  in  chemical  equilibrium. 
Accordingly  for  the  ith  component, 


or 


dui  ■  dui 


98 


Eliminating  the  chemical  potential  between  eq.  98  and  eq.  19 
results  in  the  following 1 


^i  dP!  =  ^i  dPr  99. 

x£  1  xr  1 


Since  V»/x£=MW./(  J>' Nq )  and  Vr/x£  =  MWi/(^No)  where  MV^  is  mole¬ 
cular  weight,  Nq  is  Avogadro's  number,  and  is  density,  eq.  99 
can  be  re-written  as 


Assuming  an  incompressible  liquid  (  is  constant)  and  an  ideal 
gas  so  that  P?  =  j^'RT,  the  above  equation  can  be  integrated  from 
the  reference  equilibrium  vapor  pressure  at  which 
conditions  at  the  critical  size  where  P.*=P_  and  Pi,=P!V.  The 

10  1  1C 

result  of  this  integration  is  the  followings 


pic  =  pi.  ej'P(-^'(rre-Po)/(£,PIe)) 


100 


As  in  reference  5#  "the  exponential  term  is  expanded  in  a  Taylor's 
series  and  the  first  three  terms  retained.  Eq.  100  can  then  be 
put  in  the  following  formi 


where 


Pic  =  (Pie-po>tti  +  Po 

pj  pu  2 

«i  -  a-  f  +  > 


101a 

101b. 


Eq.  101  provides  a  good  approximation  to  the  gas  pressure 
of  the  ith  component  in  the  critical  sized  bubble.  For  the 
liquids  used  here  >  .9  at  their  limits  of  superheat  (e.g.,  for 
water  at  30 5C.  «  .979).  Therefore , the  component  gas  par¬ 

tial  pressures  in  the  bubble  are  nearly  the  same  as  their 
equilibrium  vapor  pressures. 

The  equations  given  in  Appendix  II  relate  the  equilibrium 
vapor  pressure  of  the  liquids  used  here  to  temperature.  These 
equations  were  obtained  by  curve  fitting  data  available  in 
various  handbooks  of  thermodynamic  properties. 

When  data  were  available,  liquid  and  vapor  densities  of  the 
hydrocarbons  were  fit  to  a  least  squares  formula.  When  not 
available,  was  set  equal  to  .915  as  recommended  in  reference  11. 
In  either  case  the  correction  given  by  eq.  101  has  a  small 
effect  on  predicting  nucleation  temperatures  using  equations  83a 
and  90. 


3.  Nucleation  Rate  and  Description  of  Solution  Method 
Equations  83a  and  90  can  be  used  to  determine  the  nucleation 


rate  at  a  given  temperature.  This  is  a  straightforward  cal¬ 
culation  because  J  is  an  explicit  function  of  temperature. 

However,  to  calculate  a  temperature,  equations  83a  and  90  must 
be  solved  by  trial  and  error.  In  addition  it  is  necessary  to 
know  what  nucleation  rate  corresponds  to  given  experimental  con¬ 
ditions.  This  can  be  adequately  estimated  if  it  is  assumed  that 
a  droplet  will  boil  if  just  one  nucleus  of  critical  size  is 
formed  in  that  droplet. 

The  usual  procedure  (refs.  5  and  11)  is  to  assume  that  a 
collection  of  N  drops  is  heated  at  a  constant  rate  dT/dt  (deg/sec) 
This  canncontrolled  experimentally.  (Even  if  dT/dt  is  some 
function  of  time,  little  error  will  result  by  using  a  suitable 
average  value).  Since  J  is  the  number  of  nuclei  of  critical 
size  formed  per  unit  volume  and  time,  and  V  is  the  volume  per 
drop,  the  following  rate  equation  can  be  written 1 

JNY  +  dN/dt  =0  102. 

The  first  term  represents  the  number  of  exploding  drops  per  unit 
time,  and  the  second  term  represents  the  change  in  the  number  of 
drops  due  to  these  explosions.  Because  dN/dt  =  (dT/dt)dN/dT, 
eq.  102  becomes  1  *  *  -va 

JdT  =  (- (dT/dt )/V)dN/N  103 . 

The  numerical  solution  of  equations  83a  and  90  to  be  pre¬ 
sented  shortly  indicate  that  over  a  wide  range  of  nucleation  rates 
the  temperature  variation  of  J  is  approximately  exponential 
(i.e.,  a  plot  of  InJ  vs.  T  is  a  straight  line).  It  is  therefore 


a  good  approximation  to  take 


J  «  Jo  exp( (dlnj/dT)T)  104, 

where  dlnJ/dT  is  a  constant.  Substituting  equation  104  into  103 

I 

and  integrating  between  N=N0  and  T=TQ  (initial  number  of  drops 

at  T=T  )  to  N=Nq  and  T=T  (the  number  of  drops  which  have  not  yet 
o  s  s 

exploded  at  temperature  T  )  gives, 

s 

J(Tg )  -  J  (Tq  )  =-  ( ( ( dlnj/dT )  dT/dt  )/V )  In  (N/No ) . 

J  is  a  very  strong  function  of  .temperature  so  that  ^ need. be  only 
a:  few  degrees  lower  than  for  J(T_)  »  J(T  ).  Incorporating.- 
this  approximation  in  the  above  equation; gives,  r.ric  r.a;:- 

h;V  •  J  Z  m~Z  ;  '  J  c  -  /  “j  :  v"5  "  - 

J  w  -  ( ( ( dinj/ dT )  dT/ dt  )/V )  In  ( Ng/Ng )  • 

This  equation  shows  that  only  a  temperature  at  which  a 
percentage  of  drops  are  left  can  be  estimated  from  nucleation 
theory.  This  temperature  is  defined  at  the  point  where  50%  of 
the  drops  have  exploded  (refs,  5  and  23),  so  that 

J  «  .693  (( (dlnj/dT) dT/dt )/V)  105. 

The  slope  of  the  InJ  vs.  T  curve  is  determined  by  solving 
equations  83a  and  90  for  a  range  of  nucleation  rates,  and  then 
determining  the  rate  corresponding  to  the  given  experimental 
conditions  (i.e.,  V  and  dT/dt).  The  temperature  corresponding 
to  this  rate  can  then  be  estimated  from  the  tabulated  results  of  the 
numerical  solution. 

For  numerical  analysis  it  is  convenient  to  write  equations 
83a  and  90  in  the  following  formi 


-66- 


J  =  C  exp(-w/kT) 


106. 


The  values  of  C  and  W  for  equations  83a  and  90  are  the  following! 


1)  for  equation  83a, 


107a  1 


Or  128B<J£): 


_  16fl°b 
3(P"-P0)2 

Aa=  4ir  (3MW  /(4nJiN  ))2/3 

w 


107b  1 


107c 


2)  for  equation  90 » 


nk(?^b) 


l6fl'or 
W  =  D 

3(P£-Po)2 

Equation  106  can  be  solved  for  T  to  give 


k(ln(j/C)) 


108a  1 


108b. 


109. 


Expressing  P  in  atm,  J*  in  g/cm3,  m  in  g,  c  in  dyne/cm, 
and  T  in  deg  K,  equations  107  and  108  can  be  re-written  as  fol¬ 
lows  (where  N&=  jtNQ/MWa,  =  £NQ/MWb,  ma=MWa/NQ  and  mb=MWb/N0 ) » 

1)  equation  107 » 

C  *  8.786xl025((^/MWa)2/3-»-(^/MWb)2//3)(A+Bl)/(B^)^  110a 


1 i 03x10  (P"-P 


B1  = 


50,265 


110c 


1.632x10“ 


110d| 


2)  equation  108, 


C  =  3.73x10 


351  4qJb 


*2  \* 


1.632x10 


<p£“Po>' 


MW^B 


“11  ^ 


111a 


111b, 


Both  (P"-P_)  and  (PM-P  )  in  the  above  equations  are  related  to 
&  o  o 

the  equilibrium  vapor  pressure  by  eq.  101. 

The  terms  A  and  B1  in  eq.  110a  were  grouped  according  to 
eqs.  110b  and  110c  to  permit  a  comparison  of  their  relative  orders 
of  magnitude  in  the  numerical  solution.  It  may  be  noticed  that 
even  when  A  and  B1  are  the  same  order  of  magnitude,  setting 
Acl-Q  still  has  a  negligible  effect  on  the  nucleation  rate.;. 

Equation  109  with  equations  110  or  111  was  solved  by  what 
may  be  called  a  modified  "binary  chop"  method  (ref.  24).  In  this 
method  a  temperature  is  assumed  and  the  right  hand  side  of  eq.  109 
is  evaluated  at  intervals  of  IOC  until  a  change  of  sign  occurs. 

The  increment  is  then  reduced  by  an  order  of  magnitude  and  the 


right  hand  side  re-calculated  at  intervals  of-lC  until  another 


sign  change  occurs,  and  then  at  intervals  of  +  .1C,  etc.  Admittedly 
this  method  is  one  of  the  most  unsophisticated  and  slowly  con¬ 
verging  root  finding  algorithms  available,  at  least  compared 
to  the  secant  method  or  reguli-falsi  method  (the  Newton-Raphson 
method  was  ruled  out  because  the  required  derivative  was  too 
tedious  to  obtain).  It  is,  however,  easy  to  program  and  always 
converges. 

4.  Nucleation  in  Pure  Liquids  -  Results 

The  program  used  to  solve  for  the  homogeneous  limit  of 
superheat  is  shown  in  plate  1.  Equations  109  and  111  were  pro¬ 
gramed  directly  and  can  be  identified  in  the  program  listing. 

(The  particular  program  listed  is  for  a  n-decane  calculation). 

The  only  input  required  for  the  program  is  P'.'  ,  a.  .,  * ,  and  . 

Plates  2  to  7  list  the  numerical  solutions  of  the  nucleation 
temperatures  to  the  nearest  ,1C  for  nucleation  rates  ranging  irom 
J=l/cm^-sec  to  J=10^°/cm^-sec.  Listed  also  are  C  (l/cm^-sec), 

O  (dyne/cm),  W(dyne-cm),  bubble  radius  R(cm),  bubble  volume 
V(cm^),  and  number  of  molecules  in  the  nucleus  x. 

It  is  seen  that  T  varies  less  than  1C  per  two  orders  of 
magnitude  change  in  J,  and  that  there  are  between  100  and  200 
molecules  in  the  critical  sized  bubble. 

A  comparison  of  the  experimentally  measured  temperatures 
listed  in  Table  I  and  the  numerical  solutions  given  in  Plates 
2  to  7  shows  excellent  agreement  for  all  nucleation  rates  assumed 
in  the  theoretical  calculations.  It  is  , nevertheless,  useful  to 
estimate  what  rate  corresponds  to  given  experimental  conditions. 


Injected  drop  diameters  for  the  pure  fuels  ranged  from 
about  ,1mm  to  1mm  with  a  median  diameter  of  .5nim  being  repre¬ 
sentative.  The  rate  of  rise  of  these  drops  was  about  5cnv/sec 
in  the  region  of  explosive  boiling  (though  a  precise  value  is 
difficult  to  measure  ).  For  a  linear  gradient  of  roughly 
3C/cm  the  heating  rate  is  about  15C/sec.  The  value  of  dlnj/dT 
was  determined  graphically  from  a  plot  of  InJ  vs.  T  (fig.  12 
shows  such  a  plot  for  the  special  case  of  decane).  As  shown 
in  fig.  12  j  is  well  represented  as  an  exponential  function  of 
temperature  which  justifies  the  approximation  of  eq.  104. 

For  decane,  dlnj/dT  w  6.4.  Similar  plots  for  the  other  liquids 
tested  would  show  that  6<dlnJ/dT<8  which  is  consistent  with 
values  reported  in  the  literature.  An  average  value  of 
dlnj/dT  «  7  will  be  used.  From  eq.  105  the  value  of  J  which 
characterizes  the  experiments  run  here  for  homogeneous  nucleation 

c  4  6 

is  J«l,4xl0  ,  or  more  generally  10  <J<10  , 

Table  IV  lists  the  experimentally  measured  '  limits  of 

4  6 

superheat  and  the  theoretical  estimates  at  J  =10  and  10  .  Also 
included  is  the  experimental  result  from  ref.  13  for  water. 

The  large  discrepancy  between  theory  and  experiment  for  .water 
has  been  discussed  in  the  literature  (refs,  9  and  13). 

5.  Nucleation  of  Fuel  Emulsions  -  Results 

Plates  8  and  9  list  the  program  for  determining  the  nuclea¬ 
tion  temperatures  of  the  emulsions  using  equations  109  and  110. 

The  same  solution  method  as  before  was  used,  and  the  results  of 
the  theoretical  predictions  at  nucleation  rates  from  J=1  to  J=10^0 
are  listed  on  plates  10  to  14, 


For  all  nucleatior.  rates  assumed,  agreement  between  theory 


and  experiment  (Table  :i)  is  rather  good.  This  agreement  was 
not  necessarily  expected  both  because  of  the  experimental  problems 
of  measurings  emulsicn  nucleation  temperatures  and  the  assumptions 
used  in  deriving  the  theoretical  model. 

Injected  drop  i  had  to  be  greater  than  about  1mm  diameter 
to  avoid  complete  dissolution  of  the  water.  The  data  listed  in 
Table  II  were  the  ’efore  for  droplet  sizes  between  1mm  and  2mm 
diameter.  A  med  an  size  of  1.5mm  is  assumed  with  a  corresponding 
volume  of  .0141;  P,  The  rising  drops  contained  15#  by  volume  of 
water.  Though  ■*  ae  water  is  ideally  considered  to  be  dispersed  as 
tiny  droplets  ia  the  fuel,  this  is  probably  not  the  case  .  rby  the 
time  they  rea  ;h  their  nucleation  temperatures  as  suggested  in 
the  illustration  shown  in  fig,  11  .  Some  coalescence  of  the 
water  probably  occurs  so  that  by  the  time  the  injected  drop 
reaches  it.'i  nucleation  temperature,  the  internal  phase  may  be 
considered  as  one  large  body  of  water.  Whether  this  is  in  fact 
the  case  is  not  too  important  because  nucleation  at  a  liquid- 
liquid  interface  is  being  described  and  whether  this  interface 
is  betwf.  an  a  large  body  of  water  and  the  fuel  or  the  curved 
surface  of  a  tiny  dispersed  water  drop  and  the  fuel  should  not 
make  a  difference.  Hence  the  dispersed  water  volume  is  approximately 
.15x. )l4l  =  V  =  2.12xl0~3cm3. 

The  value  of  dlnj/dT  was  graphically  determined  from  the 
slope  of  a  InJvs  T  curve.  Such  curves  ave  shown  in  fig.  12 
for  all  emulsions  tested  (data  obtained  from  plates  10  to  14), 

Also  included  is  the  data  for  a  water-freon  system.  It  should 
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be  noted  that  the  experiments  reported  in  ref.  4  were  performed 
by  injecting  pure  water  drops  in  a  heating  medium  of  freon  E9» 
and  not  preparing  a  freon-water  emulsion  and  then  injecting  a 
drop  of  the  emulsion  into  a  third  fluid  as  done  here.  Their 
results  are  included  here  only  because  the  surface  and  interfacial 
tensions  of  freon  and  water  satisfy  eq.  80.  In  any  case  it  is 
seen  that  for  all  the  emulsions.; the  value  of  dlnj/dT  lies  between 
2.5  and  3.58  with  a  mean  of  about  3.  This  value  is  to  be  con¬ 
trasted  to  dlnj/dT«7  for  homogeneous  nucleation  in  the  fuel. 

The  measured  rates  of  rise  of  the  emulsified  drops  were 

around  8cm/sec  in  a  temperature  gradient  of  about  2C/cm  so  that 
dGjat  «  l6C/sec.  Therefore,  from  eq.  105 

J  -  (. 693x3x16 )/2.12xl0“3  -  1.57xl04 

or*  . 

J  *  104. 

Table  V  lists  the  esperiraental  and  theoretical  predictions 

jj, 

(for  J=10  )  of  the  nucleation  temperatures  of  the  emulsions 

tested.  The  theoretical  predictions  are  in  excellent  agreement 

with  the  measured  temperatures.  It  must  be  remembered,  though, 

that  only  the  highest  recorded  temperatures  are  listed  in  Table  V, 

The  great  majority  of  rising  drops  had  air  bubbles  attached  to 

them,  and  others  exploded  at  lower  temperatures  probably  due 

to  motes  or  imperfect  wetting  at  the  water-fuel  interface, 

A  relative  comparison  of  the  theoretical  and  measured 

temperatures  for  both  the  pure  liquids  and  emulsions  examined 

here  is  shown  graphically  in  fig.  14,  The  good  agreement  of  both 


homogeneous  nucleation  (already  well  established)  of  the  pure 

liquids  and  the  assumed  mechanism  of  droplet  explosions  in 
emulsified  fuels  is  apparent. 

The  pre-exponential  term  for  nucleation  in  fuel  emulsions 

25  1  32 

is  C  *  10  J ,  while  for  homogeneous  nucleation  C  »  1CK  .  This 

reduction  is  due  primarily  to  considering  molecules  at  the 

interface  rather  than  molecules  within  the  bulk  of  the  drop  as 

potential  nuclei  (i.e.,  N— »N2//^) 

The  values  of  A  and  B1  (defined  in  equations  110b  and  110c) 

are  both  approximately  the  same  order  of  magnitude  for  all  emulsions 

tested.  This  might  at  first  imply  that  including  the  cross 

sectional  area  of  an  a  molecule  in  the  pre-exponential  term  of 

eq.  83a  is  important.  However,  the  large  multiplicative  factor 
25 

of  10  J  in  eq.  110a  has  a  damping  effect  on  additive  terms  of 
the  same  order  of  magnitude.  If  Z)ais  set  equal  to  zero  in 

eq.  110a,  the  solution  of  eq.  109  is  not  effected  and  the  same 
temperatures  as  listed  in  Table  V  result.  This  would  lend  sup¬ 
port  to  the  use  of  eq.  71  instead  of  eq.  73  for  the  general 
expression  describing  nucleation  at  a  liquid-liquid  interface. 
However,  it  would  be  conceptually  difficult  to  explain  how 
molecules  of  either  liquid  could  penetrate  the  interface  if 
4otss  A*=  0.  This  is  so  because  the  assumption  was  made  that 
the  interface  i-s  a  boundary  across  which  molecules  cannot 
cross  unless  there  is  an  "opening"  for  their  passage.  Of  course 
this  is  not  strictly  correct.  Nevertheless,  for  purposes  of 

describing  nucleation  at  a  liquid-liquid  interface  this  assumption 
does  lead  to  results  consistent  with  experimental  measurements. 


The  vapor  pressure  of  water  is  approximately  two  orders 
of  magnitude  greater  than  the  fuel  at  the  predicted  nucleation 
temperature.  The  fundamental  condition  of  applicability  of  eq. 
83a  is  therefore  satisfied.  Even  for  a  decane-water  emulsion 
where  P"  «  28atm  and  4atm,  the  predicted  temperatures  are 
in  quite  good  agreement  with  the  highest  measured  explosion 
temperature  .  However,  just  .how  much  more  volatile  one  liquid 
must  be  than  another  for  eqs.  66  or  67  to  be  valid  is  an  open 
question  at  this  time.  The  results  presented  appear  to  be  valid 
when  there  is  at  least  a  one  order  of  magnitude  difference  in 
vapor  pressure. 

6.  Application  to  Emulsified  Fuel  Combustion 


Disruptive  burning  of  emulsified  fuels  is  well  established 
in  the  literature.  The  cause  of  this  phenomenon  is  here  attributed 
to  nucleation  at  the  internal  phase  (water)  -fuel  interface. 

Whether  droplet  boiling  will  or  will  not  occur  in  a  given  physical 
situation  is  determined  primarily  by  the  surface  tension  of 
the  fuel  and  the  vapor  pressure  of  the  water.  The  addition 
of  an  emulsifying  agent  to  the  fuel  during  emulsion  preparation 
changes  the  mechanism  of  nucleation  from  bubble  growth  between 
the  interface  -  a  lens  shaped  bubble  -  to  spherical  bubble 
growth  tangent  to  the  interface  and  surrounded  by  the  fuel  (fig. 3d) 
Because  of  this  it  is  not  necessary  to  know  a  specific  value 
of  the  interfacial  tension  of  the  two  liquids  composing  the 


emulsion.  The  following  proportionality  based  on  eq.  83a  can 
then  be  written i 


«2 
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To  determine  whether  a  particular  emulsion  will  bum 

disruptively  or  not,  it  is  necessary  to  know  only  whether  T  is 

s 

greater  or  less  than  the  boiling  point  of  the  fuel  at  the  ambient 

pressure,  PQ.  The  steady  state  temperature  of  a  burning  pure 

fuel  droplet  is  not  exactly  i'.ts  boiling  temperature  but  is  •  slightly 

less.  This  isAbecause  the  fuel  mass  fraction  at  the  droplet 

surface  is  not  exactly  1  but  is  slightly  less  due  to  the  presence 

of  the  products  of  combustion  (i.e.,  the  vapor  pressure  of  the  fuel 

is  not  its  equilibrium  vapor  pressure).  Nevertheless,  in  most 

cases  the  actual  droplet  temperature  is  usually  only  a  few 

degrees  less  than  its  saturation  temperature. 

The  heating  period  in  droplet  burning  (ref.  25)  provides 

evidence  that  disruptive  burning  or  fragmentation  of  the  primary 

fuel  droplet  will  occur  only  within  the  first  10$  to  20$  of  the 

total  burning  time.  The  main  question  is  still  whether  the 

nucleation  temperature  (Tg )  is  greater  or  less  than  the  steady 

state  droplet  burning  temperature. 

Tables  II  and  V  show  that  for  decane,  dodecane,  and  tetra- 

decane  emulsions  no  explosive  boiling  will  occur  during  combustion 

because  T  >  T.  (boiling  temperature)  for  these  emulsions, 
s  b 

For  a  hexadecane-water  emulsion  disruptive  burning  should  occur 
because  T  <  T.  (i.e.,  269.7  <  287.5). 


Nucleation  temperature  predictions  for  o  tfier  than  water-fuel 
emulsions  stabilized  by  an  emulsifying  agent  are  straightforward. 


One  only  has  to  subsitute  the  appropriate  properties  of  the 
internal  phase  and  ..'rir  fuel  into  f . ;  computer  program 

listed  on  plates  8  and  9  to  estimate  a  tea j erature .  One  need  not 
go  to  this  trouble  for  the  particular  case  of  an  emulsion  of 
methanol  and  three  of  the  four  fuels  used  here.  The  homogeneous 
limit  of  superheat  of  methanol  is  i860  (rjf.  15).  Since  a  liquid 
cannot  be  heated  to  a  temperature  higher  than  its  homogeneous 
nucleation  temperature  corresponding  to  ;he  ambient  pressure, 
internal  boiling  is  guaranteed  to  occur  for  an  emulsion  of 
.... :  i  dodecane,  tetradecane,  and  hexadeca.e.  Simply  stated,  the 
steady  state  droplet  burning  temperatui  2s  of  these  fuels  are  all 
higher  than  the  homogeneous  nucleation  temperature  of  methanol 
(for  burning  at  1  atm).  Such  an  observ  ition  leads  to  the  following 
generalization  for  the  prediction  of  c roplet  explosions  in  the 
burning  of  emulsified  fuels i  Explosiv t  burning, or  it  the  very  least 
internal  boiling, will  be  guaranteed  to  occur  if  the  steady  state 
burning  temperature  is  greater  than  the  homogeneous  nucleation 
temperature  of  the  internal  phase.  If  this  is  not  the  case,  the 
theoretical  formulations  derived  here  can  be  used  to  determine 
at  what  temperature  nucleation  should  be  expected. 

This  study  was  concerned  on?:  with  developing  a  theoretical 
model  for  disruptive  burning  in  e  nulsified  fuels  and  its  ex¬ 
perimental  verification  at  latr  ,  //hile  information  on  high  pres¬ 
sure  nucleation  can  be  obtained  'rora  the  results  presented  here, 


there  is  currently  no  experimental  evidence  which  couli  be  used 
in  support  of  any  subsequent  conclusions  (except  at  latm  as 
reported  here).  Before  this  work  is  finally  completed  an  experi¬ 
mental  study  pf  nucleation  at  high  pressures (>10atm)  will  be 
undertaken. 

AJdetailed  discussion  of  the  role  of  the  suspending  quartz 
fiber  in  inducing  explosive  burning  of  emulsified  fuels  will 
also  be  given  at  some  later  time.  Such  a  discussion  will  not 
require  theoretical  formulations  beyond  what  already  exist,  in  the 
literature. 

V.  Conclusions 

A  theoretical  and  experimental  investigation  of  the  explosive 

burning  phenomenon  in  emulsified  fuel  combustion  has  been  under¬ 
taken.  The  approach  was  that  of  using  nucleation  theory  to 

explain  the  mechanism  of  droplet  explosion  in  the  burning  of 
such  fuels.  The  theoretical  model  is  based  on  the  assumption 
that  a  spherical  bubble  grows  at  the  interface  between  the 
fuel  and  internal  dispersed  phase,  and  is  surrounded  almost  com¬ 
pletely  by  the  fuel  except  for  an  opening  large  enough  to  permit 
molecules  of  the  volatile  internal  phase  to  enter  the  bubble. 

The  bubble  grows  until  a  size,  known  as  the  critical  size,  is 
reached  at  which  the  bubble  is  in  a  metastable  state.  -Growth 
beyond  the  critical  size  occurs  by  the  addition  of  a  few  more 
molecules  and  results  in  a  sudden  and  rapid  boiling  of  the  emulsion 
drop. 

No  previous  observations  (to  the  authors  knowledge)  have 
been  reported  on  free  burning  emulsion  drops. 


The  experimental  method  reported  here  attempted  to  model  this 
configuration  by  measuring  nucleation  temperatures  of  rising 
emulsion  drops  in  a  heavier  heated  liquid.  This  technique 
enables  the  emulsified  fuel  drop  to  be  heated  to  a  temperature 
well  in  excess  of  its  boiling  point.  In  a  free  burning  emulsion 
drop,  however,  the  maximum  temperature  is  limited  by  the  boiling 
point  of  the  fuel. 

The  conclusions  of  this  investigation  are  given  below. 

1)  The  expression  for  the  nucleation  rate  at  a  liquid- 
liquid  interface  has  the  general  form, 

J  =  C  exp(-W/kT). 


When  P"»P" ,  the  following  forms  of  C  and  W  are  assumed 
d.  o 

depending  on  the  surface  and  interfacial  tensions  of  the  two 
liquids! 


la)  ffab>  V®b 
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.  2)  Setting  Aoj1  Afc=0  in  the  above  equations  has  a  negligible 


effect  on  the  nucleation  rate. 

3)  For  water-fuel  emulsions  stabilized  by  an  emulsifying 
agent,  the  interfacial  tension  will  be  very  low  and  condition 
2a  above  will!  be  satisfied.  The  nucleation  rate  is  then 
independent  of  the  interfacial  tension  and  depends  only  on  the 
pure  fuel  surface  tension  and  vapor  pressure  of  water. 

4)  Nucleation  temperatures  of  emulsified  fuels  can  be 
estimated  using  nucleation  theory.  These  temperatures  are 
well  in  excess  of  the  boiling  points  of  the  internal  phase 
but  less  than  their  homogeneous  limits  of ^superheat,  . .  .. 

5) The  frequency  of  nucleation  at  the  temperatures  listed 

in  Table  II  was  very  low  (being  approximately  1  out  of  10).  The 
great  majority  of  drops  had  air  bubbles  attached  to  them  which 
expanded  on  heating.  Dissolution  of  the  water  occured  during 
droplet  ascent  and  was  most  pronounced  for  emulsion  droplets 
less  than  about  ,5mm  diameter.  Nevertheless,  this  did  not 
effect  superheat  temperature  predictions. 

6)  The  location  of  bubble  formation  can  be  predicted  by 
determining  which  of  equations  76  to  81  is  satisfied  for  the 
particular  components  of  the  emulsion.  This  supports  the 
conclusions  of  reference  4. 

7)  Disruptive  burning  (i.e.,  droplet  explosions)  will 
occur  if  the  nucleation  temperature  as  calculated  using  eq.  83a 

is  less  than  the  steady  state  burning  temperature  of  the  emulsion. 

8)  If  the  homogeneous  limit  of  superheat  of  the  internal 
phase  in  an  emulsion  is  less  than  the  steady  state  temperature 
of  a  burning  drop,  disruptive  burning  will  occur.  This  is. sure 
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n-Decane 

142.3 

1?4.1 

228. 

3^.3 

n-Dodecane 

170.3 

216.2 

250. 

385.1 

n-Tetradecane 

198.4 

252.5 

259. 

421. 

n-Hexadecane 

226.4 

287.5 

263. 

444. 

Freon-E9 

1614.3 

290. 

228* 

363. 

Key* 

MWf -molecular  weight  of  fuel 
T^f-boiling  point  of  fuel  at  latm 
T  -  nucleation  temperature 
Tcf“critical  "temPera'ture  of  fuel 


data  from  reference  4 


Table  IV 


Homogeneous  Nucleationi 
Comparison  of  Theory  and  Experiment 


Liquid 

Tb 

T 

sm 

4 

J«10 

T 

SC  £ 

J=106 

n-Pentane 

36.1 

147. 

146.1 

146.7 

n-Hexane 

68.7 

182. 

183.4 

184.1 

Benzene 

80.1 

225.3 

228.9 

229.8 

n-Octane 

125.6 

241.7 

238.8 

239.5 

n-Decane 

174.1 

282.8 

281.9 

282.6 

Water 

100 

279. 5a 

305.4 

306.3 

Key  i 

T^-boiling  point  at  latm 

T  -measured  nucleation  temperature 
sm 

Tgc-calculated  nucleation  temperature 
using  eq.  90 


data  from  reference  13 


Table  V 


1  Nucleation  in  Emulsified  Fuels i 
Comparison  of  Theory  and  Experiment 


Emulsion 

T 

xbf 

T 

sm 

Ta 

sc 

T 

sm 

Decane-Water 

174.1 

228. 

230.8 

-2.8 

Dodecane-Water 

216.2 

250. 

252.5 

-2.5 

Tetradecane-Water 

252.2 

259. 

262.2 

-3.2 

Hexadecane-Water 

287.5 

263. 

269.7 

-6.7 

V.  ^  ■  1 

_  L 

-  lb 

Freon  E9 -Water 

290. 

228b 

231.7° 

-3.7 

Key  i 

T^-fuel  boiling  point  at  latm 

T  -measured  nucleation  temperature 
sm 

Tgc-calculated  nucleation  temperature 
using  eq,  83a 

^=10^ 

bdata  from  reference  4 


LIQUID  b 


UQOID  a. 

FI&.  I 

T^AN'S.FEt.  ARFAS  FOR.  'BUfcfcLF 

frWWTtt  AT  A 

LISlMD  -  UQimD  XNTt?FA^t 


LlftUiD  b 


T*n^P/\CE- 


b  woiecucr  / 


Fife-  2. 

SCHEMATIC  Pi A6-M/IA  OF  I NT^FACiAL  U'6-i«W 

BeTwe^M  two  immi  shelf  tiaaios 


F\G  3&.  F16.  3-  t 


P  Ife  •  3c  Fi6.  3d 


FIG .  3 

pos^likf  LOCATION)^  Cf  r  6?ounn 

(AT  A  LI&utQ -t«5UlE>  INTW/Kf 


fiD-A132  177 
UNCLASSIFIED 


NUCLEflTION  AT  fl  LIQUID-LIQUID  INTERFACE  AND  EXPLODING 
DROPS  IN  EMULSIFIED  FUEL  COMBUST IONCU)  PRINCETON  UNIV 
J  DEPT  OF  AEROSPACE  AND  MECHANICAL  SCIENCES.  . 

C  T  AVEDISIAN  JUL  76  PUAMS-1315  F/G  21/2 


2/2 


1 


Fife.  Hb 


Tot  a  l  c  ySrtev*  bo^oY\/ 
S  ub-c  bcuindary 


FIG-  £ 

OF  a  D20P  OF 
LIQUID  4  IN  A  LA ICrt 
MASS  OF  Ll^iD  b 


Fit  7a 


FIG.  7 

SCHEMATIC  iLUISTH/flV^  of  an 

EfAULSlOM 


THEWCOUPIC 


95- 


~  FI6  g 

SCHfM/TTI C  OV  EXPtW/lf  NiMH. 

APPAW^ 


Fife  II  ^ 

ENVUl<>(OfO  CXW  C>Uf*jM£  INITIAL 

QP 

A^cpnt 


Fife-  II  t 

VtoM* T  AT  SOWIr  lAt^ 
Tl/VNt 


S'OtoE  VISUAL  02&EWTTId*i£  PF  WATER-  FUEL 

EMOiyTJS  RISING  |N  6L/CERINF 


THE  CENTIMETER 
ESSER  CO  MAUI  in  Li  s  a 


Appendix  I 


'  Constants  for  Equation  93 » 
O'  =  dQ(l  -  T/Tc)n 


Liquid 

O' 

o 

n 

!c 

Reference 

Decane 

59.37 

1.38 

617.4 

15 

Hexadecane 

55.2 

1.33 

717. 

26 

Hexane 

52.4 

1.22 

507.4 

27 

Octane 

55.72 

1.30 

568.8 

15 

Pentane 

53.4 

1.23 

469.4 

26 

Tetradecane 

54.9 

1.32 

694. 

26 

Freon  E9 

37.41 

1.22 

636. 

4 

— 

- 

- 

- 

- 

Benzene i  cr=80 . 877- . 21 2T+1 . 21 836xlO-4T2 

Dodecane  i  <r=58 . 75- .  1 3354T+6 .  726xlO’-5T2 

Water i  35 . 547- . 377T+1 . 08lxl0’3T2-2 . 

197x10"6T3+1 

.  424x10’ V* 

Units j  a  -  dynes/cra 

T  -  degrees  Kelvin 


Appendix  II 

Vapor  Pressure  as  a  Function  of  Temperature 

Benzene i  P=exp(9 . 3-1 . 575x103/T-1 . 526xl06/T2+5 . 045x108/T3-6 . 295xl010/T4) 
Decane »  P=exp( 8 .130-1. 406x103/T-1 * 29x106/T2+1 . 33x108/T3) 

Dodecane i  P=exp( 7 . 035-321 . 3/T-l . 97x106/T2+2 ,129x108/T3) 

Hexadecane i  p=10^“3366,1//T+8,96/760  * 

Hexane »  P=exp (13.971-8 . 655x103/T+2 . 589x106/T2-5 . 449x108/t3+3 . 86xl010/T4 ) 
Octane  i  P=exp ( 15 . 1 6-1 . 1 6xl04/T+4. 31  x106/T2-1 . 06xl09/T3+8 . 871x101o/T4) 

Pentane i  P=exp(13. 56-7 . 223x103/T+1 . 8x106/T2-3. 36x10 8/T3+2 . Ilxl010/T4) 
Tetradecane*  P=10^“3*004* 375/T+8.629)^76o  * 

Freon  E9*  P=exp(24.76-1 .755x104+2.039x106/T2) 

Water 1  P=exp ( 1 2 . 409-4 . 62 6x10  3/T+l .121x10  5/T2 -4 . 272x107/T3 ) 

*  equation  obtained  from  CRC  Tables 


Units*  P  -  atmospheres 

T  -  degrees  Kelvin 


LEVEL  21Trif^:;=>^:r:  8415  DATE  =  761 7T  13/26/13 


DOUBLE  PRECISION  TL , T ,0 A, 03, OAB, DLA, DVA, DIB, DVB ,  PL4 , PLB , BATE, PA 


1,WA,W3,ALPA,ALPB,XA,X3,P,W, V  ,B,C,F, DL07,  T-EXP,  DSQRT,  A, B1  ,P,PA,P3, 

1,1V 

_ _  _ _ WRITE  (6,  100)  _  _  _ 

400  1  FORMAT  (25X,  »  LIMIT  jOF  SUgSaBEAT-PBCASSlf^^lil^^^ill^^^^lIl^^^ 

BATE=  1. 


2 


200 

3 


10 


T=400.15  • 

_I=  0  _ 

FORMAT  (/•  -  -  -  ;  NUCLEI Tio!Tl{VfE>^TDlOr3)^^^^^^^^^^^^^t-' 

01=59.37*  (  (1.-T/617.  4  )  *  *  1 .3  8) 

PLA=DEXP (8. 1296 50 47 5-1406. 690454/7-1290900. 971/ (?**2) ♦ 1 332072 18. / 

(T**3)  ) 

_DL  A  =.47  _  __ 

AV=6. 02 29 6 D2 3 
P=  (PLA-PATM)  *ALPA+FATM 

B=1  .-  (1  ./3 . )  *  (1  .-PATM/P) _ _ 

W=  1  . 6 3 1  9 8204  3  D - 1 1  *  (OA  **  3)  /  {  (P-PATM)  **21^# 

C=3  .7295379  43  D3 5* DSQ2?  (  (DLA*  *2)  *01/ {  ( W  A  *  *  3) *  B)  ) 

B=2.*OA/(  (P-PATM)  *1.0132506)' 

V=4. 18879*  (3**3) 

X=P*1 .0132506*7/(1^380250-1  6«$)_ 

IF (F-T?1 sl 15,10 
T=T+ 1 0. 


_ ___G 0  TO .  1 _ _  _ _ _  _ 

19  T=T+1.J  #  # 

20  T=T+. 1 

1=2  _  _ _ . _ 

IF  (F-T)  35, 35,  1  ■"■ ~  ~ 

35  WRI TE  (6,8)  TL ,  C  , OA,  PLA,R,V,W,X 

8  FORMAT ( 1  ?  =  *,F6.2,»  C  =  »,D10.3,»  01=  »,F6.3,»  PLA  =  * 

1  ,  F6 . 3,  ’  R  =  »,D10.3,/16X,»  V  =  »,Di0.3,»  W  =  ',010.3,*  X  = 

1F4.0) 

25  _  RATE=RATE*1.D2  _ _ _  __  _ 

IF  (RATE-5.  DIO)  2,2,26  ‘ ?~r -  -- , 


c.,p« 


Pure 


LIMIT  OF  SUPERHEAT-PENTANE 


NUCLEATION  SATE  = 
90  C  =  0.653D  33 

V  =  0.673D- 18 


0.1 OOD  01  »  9 

OA  =  3.512  PLA  =  14.715 

W  =  0.435D-1 1  X  =  162. 


=  0.544D-06 


HOCLEATION  RATE  = 
50  _  C  »  0.648D  33 

V  =  0.6 2 3D- 18 


0. 1 0OD  03  ^  'r  .  IM  • 

_  OA  =  3. 461_  PLA  =  14.877 

W  =  0.  407D-1  1  X  =  152. 


=  0.5  30D-06 


NUCLEATION  SATE  = 
10  C  =  0.642D  33 

?  =  0.  577D-18 


0.1 OOD  05 _ __ 

OA  =  3.411  PLA  =  15.042 

A  =  0. 381D-1 1  X  =  142. 


=  0. 5  16D-06 


NUCLEATION  RATE  = 
70  C  =  0.6  36  D  33 

V  =  0.534D-13 


0.1 OOD  07 

OA  =  3.361  PLA  =  15.207 

W  =  0.357D-11  X  =  132. 


=  0.503D-06 


NUCLEATION  FA^E  = 
40  _  C  =  0. 629D  33 

V  =  0.488D-18 


0 . 1  0 0 D  09  V" -  .ri  “ 

v,..  OA  *  3.303  -V-  PLA  =.J  5.402 

A  =  0.330D-1  1  X  =  122. 


P  =  ,0.4880-06 


NUCLEATION  RATE  = 
10  C  =  0.623D  33 

V  =  0.446  E-1 8 


0.10  0  D_1 1 _ _ 

OA  =  3.244  •  PLA  =  15.599 

V  =  0.305D-11  X  =  113. 


=  I.  0.474D-06 


LIMIT  OF  50PERH EA  I-HEXAN£_ 


NUCLEATION  RATE  = 
182.20  C  =  0.477D  33 

_ V  =  0 .  8  3  7  D- 1 9_ 

NOCLEATION  RATE  = 
182.80^  C_=  0. 474 D  33 

?  =  0..776D-18 


0.100D  01 

DA  =  3.288  PLA  »  13.132 

V 0.47 ID-11  I  »  153. 

0.100D  03 

-  DA  =  •  3.242  ..  PLA  =  13.267. 
M  =  0.442D-11  X  =  153. 


=  0.585D-06 


te>:0»570D-06 


NUCLEATION  RATE  = 
183.40  -  C  =  0.47 QD  33 

0.720D-18 

NUCL EAT ICNPA X E  = 
184.10  C  =  0.467D  33 

_  7  f  0. 659  D-1 8 

NUCLEAUCN  RATE  = 
184.  80  C  ?_  0.463D  33 
V  =  0.603D-18 

NUCLEATION  RATE  =_ 
185.50  C  =  0.459D  33 

7  =  0.551D-18 


0.100D  05 _ _ _ 

.  DA  =  3 .1  97  ^  PLA  -  =  13.  402 

»  =  0 .  4 1 4  D-1 1  ^  X  =:  j4  3^ 

0.1 00 D  07 

OA  =  3.143  PLA  =  13.562 

_  W  = _ 0.3 S 4 D-11 _ X  =  132. 

o .  i  o  o  d  09 '  r 

OA  =.-•  3.090..;VPLA  =L_13.  723 
V  =  0.355D-11  X  =  122. 

0  ._1 0  0  D _ 11 _ 

D A  = '  '  3.038  '  '7 PLA  =  13.896 
W  '  0. 329D-1 11  .T-  -  113.  . 


£.•=?*  ~=0 .5560-06 


=  0. 540D-06 


:=MA5.24Dt06 


III.':-: 


1 


PLATE  5 


LIMIT  OF  SOPSRH  EAX-OCTANE 


NOCLEATICH 
50  C  =  0, 

V  =  0. 


EATZ  =  0.100D  01 

29 1 D  33  OA  =  2.874  PLA  =  10.379 

121D-17  W  =  0.526D-11  X  =  167. 


=  0.661D-06 


NOCLEATION 
10  C  ■_  0. 
V  =  0. 


FATE  =  0 .  1  0  0  D  0  3  '  -Vil-A"'-- 

290D  33  OA  =  2.835_  PLA  .=^10.48 1; 

,1 12D-17  Si  =  0.  494D-11  X  =  156. 


=  0.645D-06 


NOCLEATION 
80  G  =  0, 

.  :  -  7  =  0, 


RATE  =  0. 100 D  05  _ 

287 D  33  OA  =  2.791-  PLA  =  10.800 
10  3D- 17  4  =  0.460D-11  _i  =145.; 


'=  “  O'.  6  27  D- 08 


NOCLEATION 
50  C  =  0. 

7=0. 


SATE  =0.1000  07 

285  D  33  OA  =  2.746  PLA  =  10.721 

948D-18  W  =  0.427D-11  X  =  134. 


=  0. 609D-06 


NOCLEATION 
20  C  =  0. 

7  =  0. 


SITE  =  0.  130 D  09  ’ :t 

28  3 D  33  3 A  =  2.702  ;  .  PLA  =  10.943 

870D-18  "  Si  =  0.397D-11  X  =  124. 


.=_  .0. 5.920- 06 


NOCLEATION 
00  C  =  0. 

7  =  0. 


PATE  =  0.  100D  J11 _ _ 

280 D  33  OA  =  2.651  "  PLA 

7680-18  W  =  0.  364D-11  X 


=  1 0 . 9 3 3  R  =  0. 5770-06 


plate _>  ... 


LIBIT  OF  SUPER HEAT-DECA  HE 


NUCLEATION  RATE  = 
,60  C  =  0. 197D  33 

_ =  7  =  0 . 1 58D-1 1 

NUCLEATION  RATE  = 

,  20 _ C  =  0.196D  33 

7=0. 147D-17 


0.100D  01 

OA  =  2.581  PLA  =  8.706 

W  =  0.565D-11  X  =  169. 


R  =  0.723D-06 


0 . 1 00 D  03 

OA  =  2.548  _  PLA  = 

V  =  0. 532D-11  X 


8.787, 

158. 


_«  ^0.706D-06 


NUCLEATION  RATE  = 
,90  C  =  0.  195D  33 

7  =  0. 136D-17 


0.100D  05_ _ 

OA  =  2.509'“:  PLA  = 

U  = -  0.495D-1 1  X 


8.383 

147.. 


NUCLEATION  RATE  = 
60  C  =  0.193D  33 

7  =  0.125D-17 


o.  i dob  07 

OA  =  2.470  PLA  - 

tf  =  0. 461D-1 1  X 


NUCLEATICN  PATE  = 
,30  __C  =  0.192D  33 

7  =  0.115D-17 


0.  100D  09  - 

OA  =  2.431  ..  PLA  • 

W  =  0.429D-11  X 


8.980 

137. 


9,0  7_7_ 
127. 


R  =~0.687D-06^ 


S  =  0.663D-06 


6  49D-06 


NUCLEATION  RATE  = 
,20  C  =  0.190D  33 

0.  103  D- 17 


0.  100  D_1 1 _ 

OA  =  2.382 PLA  ~ 

:  W  *  :  0. 392D-1 IX 


=  "  9.203  ~  R-=— 0.6  26D-W 


VIATEL. 


LIMIT  OF  SUPE3HEAT-BATEB 


NUCLEATION  FATE  = 
,70  C  =  0.151D  35 

__V  = _ 0 .  1 5 1 D-1 8 


NUCLEATION  RATE  = 
50  -^..C  j»  _  0.  14  9  B  35 
V  =  0.141D-18 


0.1 000  01 

OA  =  13.675  PLA  =  90.345 
=__  0.626D-1 .1  ,_=  X^=  1 5 9..^ 

0.  1 00D  03  '  '  '  -v  4  i 

OA  =,13.491-.---  PLA  =  91.353, 
tf  =  0. 590D-11  X  =  150. 


=  0.331  D-< 


=.  0.323D-1 


MOCLEATION  FATE  = 
,40  C  =  0.148D  35 

V  =  0.13 ID- 18 

NUCLEATION  R  ATE  = 
,30  C  =  0.146D  35 

V  =  0. 121D-18 


NUCLEATION  RATE  = 
20  _  C  =  0.144D  35 

V  =  0.11  ID-18 


0.100D  05  _ 

OA  =  13.234  '■  PLA  =  92.497 
W  =  0.55 ID-1.1  X  *  140.  ^ 

T.1  o6d"~07* 

OA  =  13.078  PLA  =  93.652 
W  = _ 0.  5 1  5D-11  __  _  X  =  130. 

"0.1  OOD  09 

:  OA  =  12.872  .  PLA  =  94.817 
Hf  =  0.481  D-1  1  X  =  121 . 


^707315D-< 


=  0.307D-I 


.299  D-( 


NUCLEATION  RATE  = 
20  C  =  0.14 3D  35 

V  =  0.102D-18 


0.1  OOD  11  ___  __ _ 

OA  =  1 2.643  PLA “  =  96. 125 
W  =  0.445D-11.1-  X  =  112.- 


=  0.290D-< 


IV  G  LEVEL  21 


MAIN 


DATE  =76177 


fLAjeS 


.14/17/22 


DQ  g  DL  E  PRECISION  TL ,  T  ,OA  ,  OB  ,  OAB  ,  DLA  ,  DV  A  ,DLB  ,  DVB  ,  PL&  ,  PLB  ,  BATE,  PATH 
1,HA,W8,ALPA,ALPS,XA,XB,R,tf,V,B,C,P,  DLOG  ,DEXP,DS3RT,A,B1  ,P,PA,PB,AJ 

1#AV 

MFITE  (6,  100) _  _ _ 

I  vV:  POP. MAT  (BOX ,  »  WAT2H-TETSADECASE 

3ATE= 1 .  » 

T=503.15 

_ 1=0 _ _ _ 1 _ 

'  FORMAT  (/'  )  ■  BUCLEATION  RATE  =  *  .  D 1 0T:3l 

-ir,  OB=54.9*  f  f  1.-T/69  4.  )  **1.32)  .  ;.V 

OA=135.547-.377235*T+.0010812*(T**2.) -2. 197D-6* (T**3 . ) ♦ 1 . 424D-9* ( 

1T**4.) 

_____  DLI =. 51 4067 ♦. 0029 2  3*T-4 .558 D-6*  (T»*2.) 

-f- >'■  DVA=D£XP  (-24.09765969*.  05538578*1-3. 1961 07 9 D-r5*  (T »* 2 

PLA=DEXP (12.40918-4626.  3"369/TM_f2  14978279/ lT^*27j~«y72^T787f  17(1*" 

1*3.)) 

_ _ PL  B  =  ( 1 .  /7  60  . )  *  ( 1 0  .  *  *  (-  3 0  0 4 .  3 7  5/_T_f  8 .  6 2  869  9)  ) 

PATM-1 . 0 

OAB=50. 066*. 0027247+PATM-. 1205* (T-298. 15) 

AV=6. 02296D23  _ _  _ _ = 

ALPA=  1  .  -  (DVA/DLA)  *.5*  (  ( DVA/DLA)  **2  . ) 

PA=  (PLA-PATM)  *ALPA+PATH 
-s/:-'  PB=  (PLB-P  ATM)  *  ALP  3  +  PA TM_ 

P=PA*  PB 

B=1.- (1./3.)*  (1.-PATM/P) 

DELA=12. 566371  *(  (3.  *WA/(12.  566371  *DLA*AV)  )  **.6667) _ 

A=1 . 026675563 D1  2* PA*  {  (P-PATS)  **2)  *D£LA/  (P*DSQRT  ( W - 
B1  *50.  26543  246  *  (03**2)  *?B/{P*DSQ3T (SB)  )  _ 

_  S=1 .631  982043D-1 1  *  (O B**3  )  /  (  (P-PATM)  **2) 

C=8 .786D25*  (A  +  B  1)  *  (1./DSQRT  ( (0B**3)  *B)  )  *  (  (DLA/MA)  **  .667+  (DLB/WB)  * 
1*.667) 

F=-S/(1.  39025D-  16*DL0G  (RATE/C)  ) _ _ 

R=2  .  *OB/  (  (P-PATM)  *  1 . 0 1  3  25D6  ) 

V=4 . 18879  *  (R** 3  )  -  ' '  \  ' 

,  _  IA= PA  *1. 01325  D5*V/(1.38025D-16*T)  7  — 

XB=PB*1.01325D6*V/(1. 380  2  5D-16*T) 

TL=T-273 . 1 5 

IF(I.EQ.I)  GO  TO  15  _ _ 

IF  (I.  E£».  2)  GO  TO  20  . V  ‘  :~r  v  7/  7  - 

IF  (F-T)  15,15,  10  ..  7  '7,,;.  '  ■'  T  -V  . 

”  GO  TO  1  . 

T=T-1 .0 

1*1  _ 

IF  (F-T)  1,1,19  .  ^  ^7  :  :  : 

T*T ♦ 1 • 

T=T+.1 

1=2 


V  _ _ ? LATE  q  •' 

j  I»  G  M^I»^^^m^j)ATg:=^76177^q>^:  ^lC/17/22 '/  3 


I P  (  F-T)  3 5 , 35 , 1  .:±T;;  5^vt  7 
35  WRITE  (6, 8)  TL,  C ,0A ,03 , PLA ,PLB ,R, V , XA ,XB , W , A, B1 

8  FORMAT  ( '  T  =  *,F6.2,'  C  =  ',m0.3,'  OA 


•  1F6.3,*  PLA  =  •  ,F6.  3,/16X,»  PL8  =  »,F6.3,»  P.  =  *#D10.3,»  V  =  • 

1 ,010.3, 1  '•“  XA  =  *,F5.0,/16X,  'XB  =  1 ,  F3.  0,  j/^_W #  D1 0.3,»ggg^l^g^. 

IF  (PATE-5.  DIO)  2,2,26  " 

26  STOP 

END 


© 


ig 


NUCLEATION  SATE  = 

0. 1 OOD  01 

228.50  C  =  0.197D  26 

OA  =  31.222 

OB  =  5.892  PLA 

=  26.  955 

PLB  =  3.494 

E  =  0.403D-06 

V  =  0.275D-13 

XA  =  107. 

XB  -  13.  V  - 

0.401D-1 1  A 

=  0.9406D  00  .  B1 

=  _0.  1609D  02 

NUCLEATION  RATE  = 
229.60  C  =  0.196D  26 

_  PLB  =  3.  570 

5fr:=---v=rr=?;'XB  =  12.- •:?  S  = 


0.1 OOD  03 

OA  =  30.966  OB  =  5.815  PLA  =  27.512 
R  =  0.  390D-06__  7  =  0.24BD-18  XAj?  93. 
0.370D-11  J_  A  *  £  0.9819D  00  31  =  0.1S68B02 


NUCLEATION  RATE  =  -. 
230.80  C  =  0.196D  26 

PLB  =  3.655 

_  XB  =11.  W  = 


OA  =  30.688  OB  =  5.731  PLA  =  23.129 

R  =  0.376D-06  V  =  0.222D-18  XA  =  89. 

0.339D-1JL  A  =  0.1029D  01  _  B1  =  0.  1525D  02 


NUCLEATION  BATE  = 
232.  10  _  _C  =  ._  0.  195B  26 
PLB  =  3.749 

XB  =10.  S  = 


OA  =  30.386  lOB  =  5.640  y.  PLAi  = :  2e.810M^ 
H  =  0.361D-06  V  =  0. 197 D- 18  XA  =  81. 

0.307D-11  A  =  0.1082D  01  B1  =  0.1479D  02 


PKirl 


29.558 


ni 


rz. 


■  u  -  i  ■  „■ » j  ■  .  *  .  I 


uahual 


SATER-DODECANE 


NUCLEATION  RATE  = 
249.60  _  C  =„  0.955D  25 
PLB  =  2.053 

\  XB  *  6.-i -  ff  = 


0.100D  01 

OA  =  26.299  _  OB  _*  7.321  _PLA^  =  39.272  _ 

R  =  0. 368D-06  7  V  -  0.208D-18  XA  *  112.- 

0.415D-1 1  A  =  .0.19210  01  r.  Bl  =  0.1025D  02 


N0CLEATION  KATE  = 

T  =  251.00  C  =  0.962D  25 

_  PLB  = _ 2.111 

V-'-^'XB  =  5.  ' 


H  = 


0.1000  03 

OA  =  25.971  CB  =  7.232  PLA  =  40.219 

fi  =  0. 355D-06 V  =  0. 187D-^1 8 _ XA_=_  103. 

"OZ  331D-11 


A  =  0. 201  90  01  '  '  .Bl  =F  0.1 004 D  02 


T  ■ 


NUCLEATION  BATE  = 
252.50  C  =  0.971D  25 

PLB  =  2. 176 

_ XB  = _ 5 .  H  = 


o.iooD  05 


OA  =  25.619  OB  =  7.138  PLA  =  41.253 

R  =  0.341D-06  V  =  0.167D-18  XA  =  93. 

0.349D-11  A  =  0.2128D  01 _ Pi  =  0.^9828D_01 


NUCLEATION  PATE  = 
254.10  C  =  0.981D  25 

PLB  =  2.246 

XB  =  4.  W  = 


:  OA  =  2 5 . 2  4  4  :^ ^0  b  ;  =  -  7.038  ;^:P LA.=l42 .Z3?  9 
R  =  0. 328D-06  V  =  0.147D-18 

0.317D-11  A  =  0. 22500  01 


XA  =  85. 

Bl  =  0.  9601  D  01 


ft 


T  =  255.90 


L  ..  n-p  — 

\  n  x  £»  — 

C  =  0.994D  25 

PLB  =  2.327 

XB  =  4.  W  = 


.  OA  =  24.822  ^  r  6.925 -.  "PLA!i=- 43: 673^^p 

fi.,=  0.31  3D-06  =  0. 129D-18  J'XA: 76.' 'LV.. 

0.  284  0-11 


A  =  0.2394D  01  Bl  =  0.9348D  01 


T  = 


NUCLEATION  BATE  =_ 
257.90  C  =  0.101D  26 

PLB  =  2.419 

XB  =  .4  H  = 


0.1 000.  11. 


OA  =  24.353  OB  6.801  PLA  =  4 STWGiMWl 
E  =  0.2980-06  V  =  0. 11 00-1 B  '  XA  -67VT"- ’ 

0.25 2D- 11  A  =  0.256 3D  01  .i  Bl 0. 9070D  01 


II 


PLATt,:Ji 


T 


HUCLEATION  FATE  =  0.100D  01 

259.20  C  =  0.599D  25  OA  =  24.048  OB  =  8.022 _  PL A  =  46.123 

PLB  =  1.271  B  =  0. 353D-06  V  =  0.133D-18  XA  =  1 13.  h 

_  ' '  XB  =  3.  »  =  0.418D-11  A  =  0.2610D  01  r  B1  =  •  0.6235D  01 


HUCLEATION  PATE  = 
260.60  C  =  0.612D  25 

PLB  =  1.315 _ _ 

r  :  r:  IB  =  3 -  »  '  = 

-  HUCLEATION  RATE  = 
262.20  C  =  0.626D  25 

PLB  =  1.367 

_ ___»  -  3 .  W  = 


NUCLEATION  RATE  = 
263.  90  _  C  =_  0.643E  25 
PLB  =  1.425 

XB  =  2.  M  = 


0. 1000  03 

OA  =  23.719  OB  =  7.931  PLA  =  47.194 

S  =  0.34  ID-06  V  =  0.156D-18  XA  =  104._ 

:  0. 386D-1 1  A  =  0.2736D  01  Bl  =  0.6153D.01 


0. 100D. 05 
OA  =  23.344 
H  =  0.  328  D- 06 

0.352D-1 1 


OB  =  7.826  PLA  =  48.440 

V  =  0.147D-18  XA  =  95. 

A  =  0.2887D  01  Bl  =  0.6059D 


01 


:  OA  =  22.945  OB  =  _  7.716  PLA  .= :  4 9 .  7 9 2z^^ 

R  =  0.314D-06  V  =  0.130D-18  XA  =  8*6. 

0.320D-11  A  =  0.3055D  01  Bl  =  0.5960D  01 


T 


NUCLEATION  PATE  =  0.1 00D  09 

265.80  C  =  0.664D  25  OA  =  22.499  03  =  7.593-  PLA 

_ ,,  PLB  f,.  1.49U-.  J  =  0. 300D-05_V  =  -  0.114D-18  _  , 

XB  a  2.  V  =  0.287D-11  A  =  0.3252D  01  Bl 


=  0. 5849 D  01 


NUCLEATION  PATE  = _ 0.100D_11 _ 

T  *  267.90  C  =  0.688D  25  OA  =  22.007  :  "OB  =  “  7.457-  PLA  =  537084?;^ff 
PLB  =  1.567  R=  0.  236D-06  V=  0.975D-19  XA  >£?  69.^ 
-:  -_L _ _i_ XB  =.  t  2  . B  :  =  :0.  2 5 5D-1 1  _  -  A  .=  0. 3483D  .01_  Bl.  =.  10.57270-  01 


1 


NUCLEATION  RATE  =  O.IOOD  01 
266.60  C  =  0.401D  25  OA  =  22.312 

PL3  =  0.690  H  =  0.34  ID-06 

'I  IB  =  2.  W  =  0.420D-11  A 


NUCLEATION  RATE  =  O.IOOD  03 
268.  10  C  =  0.4 15  D  25  OA  =  21.960 

PLB  *  0.718  3  =  0.330D-06 

y--;.  XB  =  •  1.  :  W  =  0.388D-11  '  A 


OB  =  8.604  PLA  =  51.997 

V  =  0. 167D-18  XA  =  113. 

=  0.3287D  01  B1  =  0.3439D  0 


OB  =  8.508  PLA  =  53.253 

V  =  0.150D-18  XA  =  104. 

^  0 i 345 1 D  -  023.4 11 D  70 1 


NUCLEATION  HATE  =  O.IOOD. 

269.70  C  =  0.431D  25  OA  =  21.585 

PLB  =  0.749  3  =  0.318D-06 

XB  =  1 .  W  =  0.356D-11  A 


OB  =  8.405  PLA  =  54.618 

V  =  0.135D-18  XA  =  95. 

=  0 . 363 3D  01  B 1  =  0.  3380D  01 


NUCLEATION  RATE  =  O.IOOD  07  V 
271.60  C  =  0.452D  25  OA  =  21.139 

PLB  =  0.788  R  =  0. 304D-06 

XB  =  1.  H  =  0. 322D-1 1  A 


OB  .283  iM  PLAjl* :.56.-2.7&'^5^ 

V  =  0.118D-18  XA  =  86. 

=  0 . 386 ID  01  B1  =  0.3344D  01 


58.057 


NUCLEATICN  EATS  =  0.100D  01 

230. 40  __  C  =  0.126D  25  OA  =  30.781 

h--  V™  PLB  =  0.125  P.  =  0.408D-06 

IB  *T  1.  S  =  0.385D-11  A 


OB  5.517  PLA  _=  27.922 
V  =  0.285D-13  XA  =  1 1 4 .' 

*  0.8998D  00  .  B1  —  0.  2323D  00 


NOCLEATION  PATE  =  0.1 00D  03 

231.70  C  =  0.133D  25  OA  =  30.479 

_ _ PLB  =  0.132 _ P  =  0 . 394D-06 

=_x;1.  i;.  W  =./  0.354D-1  1  A" 

t  WOCLEATION  RATE  = 0.  1 OOD . 05jjll!| 
233.10  C  =  0.141D  25  OA  =  30.153 

PLB  =  0.139  R  =  0.379D-06 

XB  =  1.  V  -  0.323D-11  A 


OB  =  5.451  PLA  =  28.599 

_V  =  0.256D-18_  XA^=  _105.__ 

=  •' 0.9463D  00  Bl  ?  0.22 98 D  00 


OB  =  5.380  PLA  =  29.343 

V  =  0. 228D-18  XA  =  95. 

s  0.9988D  00  Bl  =  0.2266D  00 


NOCLEATION  RATE  =  0.100D  07 

234 . 60_^_  C  =  0.150D  25. OA  =  29.804 

PLB  =  0.147  H  =  0.363D-06 

XB  =  1 .  W  =  0. 293D-1 1  A 


03  =  5.  304  gr  PLA^i30^l553^r 

V  =  0.201D-1 8  XA  =  86. 

*  0.1058D  01  Bl  =  0.2233D  00 


31.097 


